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Résumé. La thèse présentée traite de l'application de trois techniques méta heuristiques : 
l’optimisation par essaims de particules (PSO), l’algorithme hybride PSO avec la méthode 
de recherche par motifs : Pattern search (PS), et l’algorithme Multi-objectifs Jaya 
(MOJaya) pour la résolution des problèmes d’optimisation de la répartition et de 
l’écoulement des puissances dans les réseaux électriques. 

Pour se faire, le problème de la répartition et de l’écoulement de puissance a été 
formulé comme étant un problème d’optimisation combinatoire avec contraintes d'égalité 
et d'inégalité où une méthode de pénalité est utilisée pour manipuler la violation des 
contraintes durant la procédure d’optimisation. Deux cas d’optimisation sont considérés: 
l’optimisation mono-objective et multi-objective. 

Différents types de fonctions objectives ont été adoptées: minimisation du coût de 
combustible, amélioration du profil de tension, amélioration de la stabilité de tension, 
minimisation des pertes actives et réactives de transmission, maximisation de la sécurité du 
système. 

D'autre part, afin de montrer l’impact très déterminant de la technologie FACTS 
dans le domaine de l’écoulement optimal des puissances, deux dispositifs FACTS, les plus 
usités, sont considérés dans ce travail. A cet effet, le compensateur statique de puissance 
réactive (SVC) et le compensateur série commandé par thyristors (TCSC) sont insérés par 
le modèle à modification de la matrice d'admittance nodale et intégrés dans les équations 
de l’écoulement de puissance. 

Les méthodes proposées ont été validées sur des réseaux test IEEE-30 bus et les 
résultats obtenus sont comparés à ceux obtenus par d'autres méthodes et publiés 
récemment dans diverses références. Le traitement et l’analyse comparative des divers 
résultats montrent, en termes de qualité et de quantité, nettement et effectivement, la 


robustesse des méthodes proposées. 


Mots clé : Optimisation de l’écoulement de puissance (OPF), optimisation par essaims de 
particules (PSO), Pattern search (PS), PSO-PS, algorithme MOJaya, optimisation mono- 
objectif, optimisation multi-objectif, FACTS, SVC, TCSC, réseaux électriques. 


Abstract. This thesis represents the application of three metaheuristic techniques: particle 
swarm optimization (PSO), hybrid particle swarm optimization and pattern search (PS) 
algorithms, and À Multi-objective Jaya (MOJaya) algorithm for solving optimal power 
flow (OPF) problem in power system. 

In this work, the optimal power flow problem has been formulated as a 
combinatorial optimization problem with equality and inequality constraints or a penalty 
method is used to manipulate the constraint violation during the optimization procedure. 
Two optimization cases are considered: mono-objective and multi-objective optimization. 

Different types of objective functions have been adapted namely: minimization of 
fuel cost, minimization of voltage magnitude deviation, voltage stability enhancement, 
minimization of real and reactive power losses, and maximization of power system 
security. 

On other hand, in order to demonstrate the impact of FACTS devices in the optimal 
power flow domain, two most popular kinds of FACTS devices are considered in this 
work. Thyristor controlled series compensator (TCSC) and static var compensator (SVC) 
are modeled by the nodal admittance matrix modification model and integrated into the 
power flow equations. 

The proposed methods are carried out on the standard IEEE 30-bus system, and 
obtained results are compared to other methods reported in the literature. The comparison 


demonstrates the superiority of the proposed methods in terms of solution quality. 


Keywords: optimal power flow problem, particle swarm optimization (PSO), pattern 
search algorithms (PS), PSO-PS, Multi-objective Jaya (MOJaya) algorithm, mono- 
objective optimization, multi-objective optimization, FACTS, SVC, TCSC, power system. 
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Introduction. 


L'énergie électrique, produit devenu désormais vital, représente un enjeu 
déterminant dans toutes les formes de la dynamique pluridimensionnelle des sociétés. Cet 
enjeu pose une problématique sur, essentiellement, trois plan; les ressources énergétiques 
primaires et leur traitement, les processus de conversion et de production, les structures de 
transport et de distribution...Les lieux d’implantation des centres de charge, de fait très 
dispersés et les lieux d’installation des sources d’énergie qui doivent se conformer, en 
même temps, au site des sources énergétique primaires et à celui des centres de charge, 
impliquent des structures d’interconnexion et de distribution très déployés et très 
complexes. Cette complexité est encore plus annoncée par l’envergure du territoire et par 
la démultiplication des centres de charge. En effet, pour répondre à l’appel en énergie des 
centres de charge ainsi dispersés et encore plus diversifiés, les réseaux électriques doivent, 
dans ces conditions, se conformer à des critères pluridimensionnels ; structurels, 
fonctionnels, hiérarchiques, technico-économiques... On cite, par exemple et 
respectivement, les réseaux à structure maillée, magistrale, radiale ; les réseaux 
d’interconnexion interterritoriale, de transport territorial, d’alimentation régionale, de 
distribution locale, les réseaux hiérarchisés en niveau potentiel...En règle générale, les 
systèmes électriques sont soumis à l’effet des variations dynamiques et aléatoires des 
charges ainsi qu’à des diverses perturbations. Ainsi, le souci d’atteindre et de réaliser des 
relations structuro-fonctionnelles rationnelles justifie la recherche des états de régime 
optimaux. Les régimes optimaux des systèmes électriques sont, essentiellement, tributaires 
de la performance des éléments qui le composent ; tels, par exemple, les centrales, les 
réseaux électriques, les lignes d’interconnexion, les diverses sous stations, les divers postes 
d’interconnexion... Les éléments qui composent les systèmes électriques sont des supports 
physiques à des processus d’écoulement électro énergétiques complexes. Dans la recherche 
des états de régime optimaux, il est nécessaire de tenir compte de deux aspects essentiels 
de la problématique ; 

e Le choix rationnel adéquat des éléments qui composent le système 
électrique ; lequel choix a lieu en temps différé par rapport à leur 
fonctionnement, 

e La contrôle (gestion) optimal de leur régime ; lequel contrôle à lieu en 


temps réel par rapport à leur fonctionnement. 


Le traitement de la problématique sur le plan du premier aspect permet d’installer la 
plateforme de base pour l’optimisation des régimes de fonctionnement du système ; 
soumis, alors, à diverses contraintes et perturbations. En effet, dans une optique 
prévisionnelle, le choix d’une structure convenable, d’une configuration de ligne adéquate, 
d’un poste de transformation à nombre d’éléments optimal, d’un niveau potentiel 
conforme..…, permet de résoudre, dans une profondeur respectable et en temps différé (déjà 
pendant l’élaboration du système), une partie de la problématique de l’optimisation, posée 
dans un cadre général. Le traitement du deuxième aspect cité exigera ensuite des moyens 
complémentaires pour le suivi et le contrôle des régimes ; lesquels régimes porteront alors, 
à ce stade, un caractère complexe, aléatoire, non linéaire,... Le thème traité dans ce 
travail s’inscrit dans ce cadre global et concerne particulièrement l’optimisation des 
régimes par le contrôle de la répartition des flux de l’énergie électrique dans les systèmes. 
En fait, par le contrôle de la répartition des puissances active et réactive, respectivement, 
au niveau des centrales et sur le plan du réseau ; en utilisant les sources complémentaires 
(compensateurs FACTS) de la puissance réactive. 

Le travail réalisé est hiérarchisé en cinq (5) chapitres ; 

Dans le premier chapitre, on pose de manière consistante et très descriptive la 
problématique de l’optimisation des régimes par le contrôle de la répartition des puissances 
active et réactive dans les systèmes. On y trouve décrit le concept technico-économique 
utilisé, en règle générale, comme critère de base d’optimisation. 

Dans le deuxième chapitre on traite des moyens classiques de base utilisés pour la 
réalisation des régimes requis. Notamment, les compensateurs statiques ; leurs propriétés, 
leurs caractéristiques, leurs divers impacts... On y insère, exclusivement, la description des 
modèles des FACTS, utilisés actuellement dans les systèmes électriques. Leur utilisation 
est dictée par leurs propriétés caractéristiques et de là par les performances de leur 
contrôle. 

Dans le troisième chapitre, on présente le formalisme mathématique descriptif des 
mises en équation des diverses fonctions et systèmes à optimiser ; ainsi que les diverses 
contraintes, leur sens et les conditions de leur utilisation. 

Le quatrième chapitre décrit les méthodes et les procédures algorithmiques 
d'optimisation utilisées pour le contrôle des flux énergétiques (écoulement des puissances). 
On y trouve la proposition et l’application propre à ce travail. 

Dans le cinquième chapitre, on présente les résultats de simulations de l’approche 


et la procédure propres à ce travail. La mise en évidence des résultats correspondants par 


2 


rapport à ceux obtenus par d’autres approches (méthodes) dévoile, à travers une analyse 
comparative bien référée, le bien fondé et la robustesse de la proposition engagée. 


En conclusion, on donne le déduit essentiel du travail réalisé. 


CHAPITRE 1. PROBLEMATIQUE DE L’OPTIMISATION 
PAR LES FLUX DE PUISSANCE. 


1.1. Introduction. 


Les états de régime des systèmes électriques sont fondamentalement tributaires de la 
variabilité de leurs charges et de leur déploiement. Cette variabilité implique celle des flux 
de puissance de ligne ; lesquels, par leur impact déterminant justifient fortement la 
recherche des états optimaux de distribution. Les orientations dans ce sens portent, en fait 
et essentiellement, sur l’optimisation des flux de puissance, respectivement, active et 
réactive ; dont la nature physique et les diverses propriétés déterminent les lieux respectifs 
de l’application. 

Le traitement de ces processus dépend de la complexité du système considéré, aussi 
bien dans la structure et le déploiement que dans la hiérarchie et la fonction. Tel ou tel 
objectif à atteindre est soumis à telle ou telle condition (contrainte). Ce domaine est 
largement investi par des recherches diverses ; puisque la portée de son impact sur les 


aspects économiques est très importante. 
1.2. Contrôle par le flux de la puissance active. 
1.2.1. Répartition entre centrales. 


Dans la problématique globale de la recherche des régimes optimaux des systèmes 
d’énergie électrique et de leur performance, l’optimisation par la puissance active utile, à 
transmette aux consommateurs, s’avère un élément important dans les solutions requises. 
L'objectif, dans ce cas, est la répartition optimale de la puissance active générée sommaire 
sur les centrales de production sous le critère du minimum de la fonction coût de leur 
énergie primaire (coût de production de l’énergie électrique). En fait, la question est 
soumise à un traitement technico-économique au bout duquel on doit atteindre la solution 
la plus rationnelle. Ce concept est très connu, largement appliqué, puisque de base, donc 
général et pluri paramétrique. 

Les méthodes de base de traitement analytique des ces questions sont diverses 


(multiplicateur indéterminé de Lagrange, méthode du gradient...) [1, 2, 3, 4, 5, 6]. 


On reprend ici quelques considérations descriptives sous le concept de Lagrange. 
La fonction coût, dans ce cas, est introduite dans l’équilibre de la fonction de 


Lagrange, sous forme suivante, 
F=D+20 (1.1) 


où ; F - Fonction de Lagrange, D — Fonction coût de l’énergie primaire du système, À — Facteur indéterminé 


de Lagrange, 0 — Contrainte d’optimisation. 


Dans l’expression (1.1), la grandeur D représente la somme des fonctions coût de 
l’énergie primaire de toutes les centrales considérées (soit m). Ces fonctions dépendent 


évidemment de la puissance active générée (PGz) correspondante ; 
D=Y D,(P;) (12) 
k=1 


En règle générale, le coût D; est introduit par le coût spécifique d’une tonne (7) 


de l’énergie primaire (EP) et le débit horaire de cette dernière, tel que, 


D= 2 DiFa)= LeB, (2% 


où ; ck — Prix spécifique d’une tonne de EP, Bk — Débit de EP par heure, B; = f(Pcy). 


La fonction © représente la condition d’optimisation ; laquelle condition stipule 
que toutes les variables P& doivent satisfaire l’équilibre du régime permanent par la 


balance des puissances actives et s’exprime, 


O= D Por D Peur AB (1.3) 


Où ; Pcnk — Puissance active de charge dans le nœud «k» du système, 1P; — Pertes actives sommaires dans 


le système. 


La fonction de Lagrange prend, alors, la forme 


5 


m m+l 


F=È DAS PT Peau AR.) (14) 


k=1 k=1 


L’optimisation dans ce cas consiste à déterminer les variables PGz; et À pour 
lesquelles les dérivées partielles de la fonction «F» s’annulent. La recherche des 
solutions prend, ensuite, pour des spécifiques complémentaires différentes et de divers 
objectifs à atteindre, des formes différentes. Si on suppose que les pertes actives dans le 
système et la puissance sommaires active des charges ne dépendent pas de la variation des 
puissances générées (c.-à.-d constantes), alors les dérivées partielles de «F>» posées 


nulles, donnent 


= = ,....+ = —À = Cte (1.5) 


Cette dernière expression, adjointe à (1.3), défini, avec, un système de « m+1 » équations 
à «m+1 » inconnues (Pc, À ; k = 1=m). 

La dérivée partielle 0D4/0Pcx = &x, dite accroissement relatif du coût de l’EP, représente 
une caractéristique spécifique importante des centrales. Elle est fonction de la puissance 
générée et diffère pour différents types de centrale. 


La recherche expérimentale dans ce domaine à donné lieu à la forme 


2 
Ex = dog + hrEcr + Arc (1.6) 


où, les facteurs dépendent du type de la centrale considérée. 

Les pertes de la puissance actives dans le système varient en fonction de 
l’écoulement des puissances de ligne ; lequel écoulement est fondamentalement déterminé 
par la variation des puissances de charge. Si, pour l’intérêt d’une analyse quelconque, on 
suppose que la charge des nœuds du système est invariable ; alors les pertes actives ne 
dépendront plus que de la répartition des puissances actives générées des centrales. En fait, 
elles ne dépendent que des puissances de «m-l» centrales ; puisqu’une des centrales 


représente un nœud de base ; où la puissance générée est la somme algébrique de toutes les 


autres. C'est-à-dire, toutes les variations des puissances des «m-l» centrales sont 


contenues algébriquement dans celle de base (P» = PGm), tel que ; 


Ces considérations permettent de déduire, à partir de la fonction de Lagrange, 


Phases ass 
OPéz OP: 

ete 

‘P, 


(1.7) 


où, la dérivée partielle ©AP5 /0Pcx est dite accroissement relatif des pertes, suite à la 
variation de Paz. 


En posant, 


ëD, CAP. D, 
=£; —*=0,; ë, 
F4 OFGx 0, 


On réécrit, à partir des conditions d’optimalité, 


2 _3=Cte (k=1:m-D) 
A (1.8) 


&,=—-À=cte (b=m) 


Ces dernières équations, adjointes à la condition (1.3), donnent un système de « 
m+l»> équations non linéaires à « 2m » inconnues : « m » puissances Pcx (k=1=+m) ; un 
facteur de Lagrange À et «m-1 » valeurs ox (k=1-m-1). Pour résoudre ce système, il est 
nécessaire d’exprimer les pertes actives en fonction des puissances générées PG4. Une des 
procédures largement utilisée de résolution du système aux dérivées partielles (1.8), (1.3) 
consiste à introduire o4 (k=1-m-1) comme fonction linéaire des puissances générées. 


Dans ce cas, le système en question se réduit à celui à « m+1 » variables. 


1.2.2. Répartition dans le réseau maillé. 


L’optimisation par les flux de la puissance active peut être posée pas uniquement 
dans leur répartition entre centrales ; mais également au niveau des réseaux électriques 
avals. Supposons qu’il s’agisse de minimiser les pertes actives dans une ligne (L,) radiale, 


(Fig.1.1). 


2 Fig.1.1. Schéma de principe. 


52 


Les pertes actives de puissance s’expriment, 


2 2 

_ +Q 
= 3 fr 
U, 
Pour raison de maintien de la tension dans les limites admissibles, on installe un 
compensateur de puissance jQc au niveau du nœud « 4 » de charge. Les pertes actives se 


ré expriment, alors, 


apr _ À (Qi -0c) 


V 12 v° 
U, 


La variation de Qc permet d’atteindre des valeurs de AP : plus ou moins grandes, 


voire minimales quand Qc = Q4. 


ap À 


T 
V.min 2 V 
U, 


Dans ce cas le facteur de puissance est égal à l’unité cospy =1 (compensation 
totale). La minimisation des pertes actives exige un investissement en moyen de 
compensation qui se traduit par des dépenses. En règle générale, les régimes les plus 
rationnels des systèmes électriques sont ceux qui obéissent à la minimisation des dépenses 
rapportées annuelles ; ces dernières sont traitées conformément à des critères technico- 
économiques déterminés. Le cas simple de la ligne Ly radiale considérée ne contient pas 
l’effet mutuel d’autres lignes, comme c’est le cas des réseaux électriques maillés, où le 
problème de la minimisation devient plus complexe. 

Considérons, pour l’exemple, le cas de la ligne Lry dans la boucle, (Fig.1.2). Le 
courant (puissance) de la ligne ZLyy dépend des paramètres de la boucle. Autrement dit, 
toute variation des paramètres de la boucle implique une nouvelle répartition des 


puissances de ligne, donc de pertes actives. 


Fig.1.2. Répartition à minimum de pertes 


Les pertes actives sommaires dans une boucle fermée s’expriment 


où ; Px O4 -— Puissance, respectivement, active et réactive de ligne, k — Numéro d la branche, r,; — Résistance 


active de la branche, U, — Tension nominale du schéma rapporté (réduit) à un niveau choisi. 
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Il s’agit de déterminer les puissances, respectivement, active P, et réactive Q7 


correspondantes au minimum des pertes actives sommaires dans la boucle ; c’est-à-dire de 


définir un état de distribution optimale des puissances de ligne du point de vue pertes 
actives. 


La circulation des flux de puissance de ligne provoque les pertes, 


AP, =] (PF +07 )r: 4(PÈ +05) +( PE + O5 )rin + (PE + O8 Jr ]. 
N 


La dérivée par rapport à P, donne, pour le régime extrémal, 


Pr = EP, 
Py =E+P=E+E -P 
Py =P+Py=E+E+E-EP 


On réécrit, donc, 


2Pr, -2(6 Pr -2B +, —Pry 2 +8, +8 -P)ry =0 


ce qui donne, 


BG + tn +) EC + 


Pin + lv) = Bin + y) By = 0 
ou bien, 
Po = Bio — By — Ey = 0 
donc, 


P 2 Ero + Bo + Bo 
T(AP=min) — 


Too 
Analogiquement pour la puissance réactive, 
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2 Qo + Bo + Do 


Oéap-nin u 
oo 
0 rI 1 Fi 2 ru 3 r] 0 
— Æ— — Æ— 
ST4Pmin) S(4Pmin) SH1(4Pmin) Sy APmin) 
S S2 S; 


30’ 


"20° 


Too” 


Fig.1.3. Puissances optimales de ligne 


et, respectivement, pour une boucle quelconque à «n » charges avec Ux =U», (Fig.1.4) ; 


n n 


Dhr Dre 
1 1 


r N 
AB 
D 
I 


Prcap=min) = 


= 
AB 
Dr 


ee ne += 


Fig.1.4. Cas général. 


1.3. Contrôle par le flux de la puissance réactive. 


Chaque régime requis optimal, en fonction des diverses variations, exige une 
répartition de puissance réactive de ligne déterminée ; c.-à.-d, exige une balance de 
puissance réactive conforme. En règle générale, cette balance diffère de celle exigée par les 


consommateurs. De ce fait, le différend entre la puissance réactive de charge et celle de 
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ligne aboutissant au nœud correspondant est nécessairement comblé par l’installation de 
sources complémentaires de puissance réactive (compensateurs). 

Le traitement des régimes optimaux par le flux de puissance réactive a plusieurs 
orientations. Dans les réseaux de distribution, on recherche plutôt à améliorer certains 
indices de qualité de l’énergie électrique, tel que «coso », pertes de tension, pertes de 
puissance ; dans ceux d’alimentation, on s’applique plutôt à la régulation de la tension, à 
minimiser les pertes. Dans les réseaux de transport, on s’oriente à atteindre (réaliser) une 
balance énergétique conforme à un plan de tension requis. Dans cette diversité d’objectifs, 
le critère de base reste, dans sa forme générale et essentielle, 1 »adéquation entre les 
quantités de puissance (énergie) réactive de contrôle à installer et leurs lieux d’installation 
conformes au coût minimal. 

Deux cas peuvent se présenter; 

e La valeur sommaire nécessaire de la puissance réactive de contrôle est connue ; il 
s’agit alors de déterminer les quantités à distribuer dans tel ou tel nœud 
conformément à un critère donné, 

e La valeur sommaire de la puissance réactive nécessaire est inconnue ; il s’agit, dans 
ce cas, de rechercher une répartition optimale conforme également à un critère 


défini ; la valeur sommaire introduite définit, ensuite, la balance énergétique. 


La structure maillée et déployée des systèmes électriques, les processus 
énergétiques non linéaires, dont ils sont le siège, la nature aléatoire des charges et leurs 
diverses caractéristiques,..., font que les procédures de traitement des questions de 
répartition optimale sont complexes. Mais, pour seulement poser la problématique, il nous 
paraît suffisant de rapporter un cas simple. 

La figure 1.5 représente un réseau constitué par trois lignes radiales alimentées à 
partir d’une sous station S/S.A de tension déterminée. Il s’agit de déterminer la 
distribution optimale des puissances réactives sur les compensateurs à installer aux bornes 
des trois charges, sachant que la puissance sommaire de compensation, définie par la 
balance requise du système d’ensemble, est constante (Q; = Cte). 


Pour ce problème, on pose donc 
es = Qu + Qu + Qu = > A 
1 
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Le problème est posé ainsi : Quelle est la répartition de la puissance réactive sur les 
sources installées aux nœuds de charge à laquelle répond le minimum des dépenses 
rapportées ? 

Toute répartition de puissance réactive sur les sources « Ok; », implique des flux 
d’écoulement correspondants différents dans les lignes du réseau. Mais, à chaque flux 
d’écoulement différent correspondent des pertes actives différentes. 

Les dépenses réduites impliquées par le seul effet de la répartition du réactif de 


compensation, s’expriment, 


D,(Q,)= D,(Q,)+D,(Q,) 


où ; la composante « D{Q%) » représente le coût (investissement) des équipements relatif à 
l’ensemble des compensateurs. Puisque la puissance sommaire de compensation installée 
est constante, alors cette composante ne dépend plus de sa répartition ; sa variation est 
nulle (sa dérivée est nulle). On peut donc ne pas en tenir compte dans le traitement 


considéré. 


Fig.1.5. Répartition du réactif à minimum de pertes. 


La composante « D,(Q;) » dépend de la répartition de Q;- sur les compensateurs ; 
puisque les pertes actives dans le réseau dépendent de cette dernière. 


Ainsi, pour le réseau donné, on écrit ; 
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2 2 < 
D, -(@-9) rrp+(@T Qu) PE 
7 < U, US, 


où ; r; — Résistance active des ligne, 7; — Temps de pertes maximales, relatif au temps d’utilisation de la 


charge maximale, b — Prix du KW.h. 


En introduisant le multiplicateur indéterminé de Lagrange, on écrit ; 


V=D. +20 
avec, 
0=Q5-(Qi + +0) (1.9) 
Ensuite, 
so ONE A 
09, Us 00, 


Ces dernières égalités permettent d’obtenir ; 


(a 0) — (9, -0,,)nt 
(a, Qi) a (Q, —Q)nts 


D'où, on peut tirer ; 


nt 
2-0 (2-01) 
Ti 
Qu =Q, (Q, Qu), 


L'égalité (1.9), donne, ensuite ; 


Ou = Our + + jre+o-onn[-+) 


Mn ED Bt DT Bt 


ce qui permet d’obtenir ; 
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et, analogiquement ; 


1 1 
Os +OnT) + Q, -Q, 
0e FT, T3 
k2 
1 1 1 
FT; + + 
Mn ED 8 
1 
Os +Qn7: + ee Q, -Q, 
1% TD 
Qu = 
1 1 1 


1.4. Le concept technico-économique. 


En électro énergétique, le concept technico-économique représente un moyen 
procédural fondamental pour le traitement des projets d’élaboration et de réalisation des 
systèmes d’énergie électrique et pour la recherche des états optimaux de leurs régimes. En 
utilise, en fait, le concept de fonction « dépenses rapportées » (fonction coût) à un temps 
unitaire de réalisation. Cette fonction intègre tous les éléments et les paramètres descriptifs 
des divers aspects technique et économique des projets (investissement, pertes d’énergie, 
type du système, niveau de tension, investissement complémentaire, variation des 
graphiques de charge, entretien, personnel...), porte, de ce fait, un caractère général et est 
pluridimensionnelle. 

Généralement, pour simplifier le traitement, on considère que l’investissement 
initial est fait en un temps déterminé, généralement une année, représentant un temps de 
réalisation et on suppose que les diverses dépenses après réalisation restent constantes pour 
chaque année d’exploitation. 


L’expression de la fonction coût, utilisée dans ce cas, a la forme générale suivante ; 
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I 
D, = + Da + Da = Palo + Dexp + Da (1.10) 


a 


où ; /, — Investissement initial, T, — Temps d’amortissement de l’investissement de départ, p, = 1/T, — 
Coefficient normatif, D. — Dépenses d’exploitation annuelles, D; — Dépenses dues à divers dommage 


causés par interruption d’alimentation. 


Les dépenses d’exploitation peuvent être représentées par deux composantes, 
D = D, +D, (1.11) 


Où ; Diny — Dépenses d’exploitation liées à divers renouvellement d’équipement, à divers entretins, au 
traitement du personnel, au transport... ; estimées comme pondération de l’investissement initial, D, — 


Dépenses d’exploitation liées au coût de l’énergie perdue annuellement. 


En tenant compte de (1.11), (1.10) se réécrit, 
D,=D,+D,+D, (1.12) 


D; - Dépenses sommaires relatives à toute forme d’investissement. 


L’expérience dans la réalisation des projets des systèmes électriques et dans leur 
exploitation montre que les solutions optimales requises ne peuvent être atteintes par, 
simplement, le traitement de (1.12) ; puisqu’on n’y tient pas compte de certains facteurs 
déterminants impliqués par la variation en fonction du temps des graphiques de charge, qui 
s’expriment par divers apports complémentaires en éléments du système, en niveau de 
tension… 

Pour rendre ces traitements plus objectifs, les projets de réalisation des systèmes 
électriques sont revus dans le concept de distribution de l’acte de réalisation sur plusieurs 
étapes. En règle générale, chaque étape est caractérisée par son investissement propre et 
implique ses propres dépenses. Le concept d’investissement progressif sur des étapes 
successives est très justifié par celui de la planification générale de l’économie nationale 


qui détermine pour chaque année la limite de son du financier. 
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ES 


Les investissements et dépenses à consentir sur les différentes étapes pour la 
réalisation des projets peuvent être déterminés à partir (par l’évaluation) de l’impact 
nuisant des dommages économiques relatifs impliqués en cas de leur manquement, [5, 6]. 

Chaque manquement à un investissement J ; c.-à.-d, un investissement n’ayant par été 


utilisé une année, implique un impact à deux composantes ; 


e Une composante D mesurant le dommage propre de l’investissement raté (non 
utilisé) dite « dommage direct », 

e Une composante D' mesurant le dommage (mutuel) induit par l’investissement 
raté dans les autres domaines qui auraient pu bénéficier de son utilisation ; dite 


« dommage indirect ». 


Ainsi, la composante directe D étant mesuré par l’investissement 7 et la 
composante D'°- par une pondération E,; sur 1, on peut écrire, pour l’année numéro 


«1»; 


D,=I+E,l=I(1+E,) 
où ; E,y — Coefficient normatif d’efficacité de l’investissement de l’année « une » de la période considérée. 
Si ce même investissement n’est pas utilisé deux ans, alors le dommage impliqué 
sera constitué par celui sommaire de la première année et celui indirect de la deuxième 


année ; 


D, =I(1+E,)+E,1(1+E,) 


En supposant que E,7 = E,y2 = … = E,, on réécrit ; 


Aüinsi peut-on généraliser si l’investissement / n’a pas été utilisé toute la période 


«T »> de réalisation ; 


D,=I(1+E,) ARE) 
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Cette égalité exprime, en fait, la nécessité du concept de la distribution de 
l’investissement dans le temps. 


On écrit, donc analogiquement, respectivement, pour les années « 2 », «3 »,...; 


D;.= 
T2 
D,,=1(1+E,) 
où ; Z1+E,)' — Facteur de rapport des dépenses sommaires (dommages) à l’année de référence. 


Pour atteindre une forme généralisée des dépenses rapportées (fonction coût ou 
fonction objective), correspondante à un manquement étendu à toutes les années de 
réalisation, on procède dans [1] par raisonnement plus simple. 

En effet, pour des dépenses réalisées sur t- années à partir de l’année initiale de 
référence, on exprime, conformément à la formule des pourcentages complexes (1.13), les 


dépenses rapportées comme ; 


D se D 


d (1+ p,) 


où ; 1/(1+p,) = p — Coefficient de rapport des dépenses à l’année de référence, t — Temps de réalisation des 
P P pp p 


dépenses. 


Pour la commodité du raisonnement, supposons que les étapes sont égales chacune 
à une année (f=]) et que l’investissement se fait au début de l’année et le compte à sa fin. 
S1 on suppose également que les dépenses de la première année restent constantes pendant 
toutes les autres, alors, les dépenses sommaires sur un nombre infini d’étapes 


d’exploitation, réduites à la fin de l’année initiale seront ; 


D D D 
D = — À +2 ++ © — 2 pD, + p'D, + p'D, +. (1.14), 
(1+p,) (1+p,) (1+p,) 
ou, alors, en supposant que D; = D: = D3 =... =D. ; 
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D, = pD,+p°D,+p°D,+..= pDi(1+p+p°+..)= er = dpD, 


où ; p = 1/(1+p,) — Coefficient de rapport des dépenses pour une année. 

Dans les projets réels, les dépenses au cours de leur exploitation ne restent pas 
constantes ; mais augmentent, en fait, en fonction de la charge et du déploiement du 
réseau ; ce dont on doit, nécessairement, tenir compte. En effet, si à partir de la deuxième 


année, les dépenses augmentent d’un écart AD, (par rapport à D) et restent, par exemple, 


constantes le long des années suivantes, c.-à.-d. ; 


D, = D, +AD, 
D, = D,+AD, = D, +AD, =D, 
D, = D, +AD, = D, +AD, = D... 


alors, les dépenses rapportées sommaires augmentent par rapport à (1.14) de 
D;;=p"D; (1+ p+p” +.) = dp’AD, 
De même, s’il y a un écart de variation à partir de la troisième année ; 
D,, = dp,AD, 


et ainsi de suite sur le délai infini d’exploitation ; ce qui permet, donc, d’écrire pour ce 


dernier ; 


où ; AD, — Différence (écart) entre dépenses, respectivement, relatives deux années qui se 
suivent, t et t-1 ; c’est une fonction pluri paramétrique. 


Si la période d’exploitation est limitée au temps T, alors ; 
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1.5. Graphique de charge. 


Les graphiques de charge représentent une donnée fondamentale pour le traitement 
des problèmes de pronostique, de planification et de diverses questions d’optimisation des 
systèmes d’énergie électrique. Relativement aux différents cas à traiter et à leurs 
caractéristiques spécifiques, on en distingue les graphiques de charge, respectivement, 
journaliers, hebdomadaires, mensuels, trimestriels, semestriels et annuels. En règle 
générale, ces graphiques portent un caractère complexe (aléatoires, non linéaires...) ; ce 
qui rend difficile leur utilisation immédiate dans les traitements analytiques. Pour 
contourner ces difficultés, on recoure souvent à des ajustements adéquats à rendre plus 
simple leur introduction, en tant que donnée, pour diverses descriptions et quantifications 
des régimes. Ces graphiques sont, en fait, habituellement définis par leur valeur, 


respectivement, maximale (valeur de pic), minimale (valeur de base), moyenne (Sy), 


(Fig. 1.6). 


S(t) 


| 


T tt 
Fig.1.6.Caractéristiques descriptives 
du graphique de charge 


Pour les analyses descriptives et estimatives, on utilise également les grandeurs comme ; 


e Le temps d’utilisation de la charge annuelle maximale, 
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e Le temps de pertes maximales annuelles, 


fs ?(t)dt 


6? 


max 


T 


e Le rapport, 


T 
[ S(t)dt 
0 

Sn se T 1 


Se Le à 
rex Î S(t)dt 
0 


T 


max 


e [La relation, 


Fi) 


Où ; À = [Stat - Energie annuelle, définie par le graphique de charge correspondant. 


Par la valeur de TX, on peut distinguer trois (3) groupes de systèmes électriques ; 


e A faible densité de charge, Tnax < 4000 h, 
e A moyenne densité de charge, 4000h< Tax < 6000 h, 
e A forte densité de charge, Tax > 6000 h 


Ces subdivisions permettent ensuite d’ajuster les graphiques typiques adéquats. 
Dans le traitement des graphiques de charge, on doit nécessairement distinguer les 
puissances produites (consommées) de celles installées. 

Les processus interactifs dans les systèmes d’énergie électrique d’ensemble 


s’expriment par la balance énergétique entre les centres de production et ceux de 
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consommation, moyennant les pertes d’énergie dans les réseaux de transport et de 
distribution. 


On en distingue, ainsi, la balance de la puissance active ; 
Dri=Zrs+XAr 
i i i 
et celle de la puissance réactive ; 


De +» (LE 1 0; => Qc; +) A0, 


Où ; 2PG; ZOc;: — Puissance sommaire, respectivement, active et réactive des générateurs ; ZPci, ZQci — 
Celles correspondantes des consommateurs ; ZAP;, ZAOQ; — Respectivement, pertes sommaires actives et 
réactives de puissance dans les réseaux ; ZPmi ZO;— Puissances réactives sommaires des compensateurs 


et des lignes. 


La différence dans le sens physique entre les ces deux grandeurs implique une 
différence d’impact de leur écoulement dans le système électrique ; autrement dit, dans la 
stratégie de contrôle des régimes par leur répartition respective. Par exemple, la balance de 
la puissance réactive diffère de celle de la puissance active, essentiellement, par le fait 
qu’un surplus de puissance réactive dans un quelconque rayon (une partie) du système ne 
peut toujours compenser un déficit dans un autre rayon. Ceci s’explique par le fait que le 
transfert de la puissance réactive sur des grandes distances n’est pas toujours rationnel et 
n’est pas toujours possible ; puisqu'il implique des variations de la tension aux bornes du 
système récepteur. 

Pour le contrôle des flux de puissance réactive, 1l est nécessaire donc de prévoir des 
réserves correspondantes ; lesquelles réserves sont représentées par des sources 
complémentaires, installées principalement dans les nœuds de charge exigeant un haut 
niveau de qualité de l’énergie électrique. Ces sources sont aussi prévues pour relever le 
facteur d puissance des turbogénérateurs et pour réduire les pertes de puissance active et 
réactive des lignes et des transformateurs. 

En règle générale, chaque « kw.h » nouvellement installé exige «0,9-+7,2 kvar » 
(AQcomp = 90 + 120 %P) de puissance réactive de compensation (de source 


complémentaire). Ce rapport d’implication entre puissances est très significatif quand il 
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s’agit de traiter la problématique du transport de l’énergie électrique sur des grandes 
portées. En effet, les réseaux ou les lignes de transport de l’énergie électrique et 
d’interconnexion régionale et interterritoriale entre de grands systèmes électriques sont 
caractérisés par des flux de grande puissance réactive, conduits sous très haute tension ; la 
puissance réactive de compensation est ainsi considérable. L’aspect, ainsi déterminant de 
ces relations dans les lignes de grandes longueur, pose un sérieux problème quant au 
positionnement optimal dans la gestion de leur régime. 

Cet aspect de la problématique posée peut faire l’objet d’une perspective et d’un 
approfondissement, orientés vers la recherche des conditions à réaliser pour atteindre une 
meilleure performance des liaisons d’interconnexion de grande portée électrique. 
Précisément, les conditions à installer doivent s’exprimer par des relations optimales 
requises entre les paramètres de configuration de la ligne et les paramètres du régime, 
conformément aux indices (critères) requis des états de régime ; tel que plan de tension, 
capacité de transport, stabilité, rendement, réserve de puissance. 

Dans tout problème d’optimisation, les systèmes à traiter sont introduits par des 
modèles à divers paramètres d’entrée et divers paramètres de sortie telle une boite noire, 
(Fig.1.7). Conformément à la nature du système (objet) traité, aux conditions de 
fonctionnement, aux critères d’optimisation.…., les paramètres d’entrée et de sortie peuvent 
être de différents types : contrôlables, non contrôlables, mesurables, non mesurables. La 
recherche de la solution à partir des équations de liaison est ensuite soumise à des 
conditions d’optimisation, introduites par des limites de type « égalité » (balance) ou de 
type « inégalité ». 

Sur la figure 1.7, on donne un exemple de modèle de principe d’un système 
d’énergie électrique constitué par « n » centrales C; électriques où, B; — représente le débit 
par centrale de l’énergie primaire, P; — la puissance active de chaque centrale, P>; — 


Puissance active sommaire de sortie du système d’ensemble. 


b. 
Pi . , ot Hu es 
Ca Q—— Fig.1.7. Schéma de principe de modélisation 
d’un système d'énergie électrique ; a. — Système 
a. d'ensemble, b. — Représentation d’une centrale. 
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En règle générale, le système à traiter peut renfermer des centrales thermiques 
uniquement ou un certain nombre de centrales de sources énergétiques primaires 
différentes. Dans ce cas, le mode de traitement de l’optimisation peut différer en fonction 
de la nature de la ressource primaire. 

Les centrales thermiques, fonctionnant dans des conditions normales, n’éprouvent 
pas de limites en réserve d’énergie primaire. Si la centrale dispose d’une réserve requise 
normalisée en carburant, elle peut fonctionner à pleine charge. Pour ce système, la période 
d’optimisation peut correspondre à n’importe quel intervalle de temps. 

Les systèmes d’énergie mixtes, contenant les centrales hydrauliques, sont 
caractérisés par e fait que ces dernières éprouvent une limite quant à l’utilisation de la 
réserve en énergie des réservoirs d’eau. L'utilisation de 1 ressource hydraulique est limitée 
au courant de la journée, de la semaine, du mois et de l’année en fonction de la période de 
réglage du stock de la centrale. De ce fait, 1l est exigé de tenir compte de la limite du stock 


sur la période du traitement. 


24 


CHAPITRE D: LES DISPOSITIFS FACTS. 


2.1. Définition et concept FACTS. 


Les FACTS (Flexible Altenating Control Transmission Systems) sont des systèmes 
de transmission du courant alternatif incorporant des contrôleurs à base d’électronique de 
puissance et d’autres contrôleurs statiques à des fins de contrôle des réseaux électriques en 
améliorant la capacité de transport des lignes, la stabilité du plan du potentiel, en 
optimisant ses régimes... [11]. 

Le concept FACTS a été introduit dans les systèmes électriques vers les années 
1988. Actuellement, ils sont devenus un outil incontournable pour l’amélioration de 
l’exploitation des systèmes et pour la flexibilité de leur régime ; laquelle flexibilité est une 
propriété importante pour la maîtrise de leur contrôle à même d’atteindre une bonne 
performance. 

Ils sont caractérisés par une commande rapide et une durée de vie respectable. Avec 
les avancées récentes dans la technologie des semi-conducteurs, le temps de réponse de ces 
dispositifs a diminué jusqu’à quelques millisecondes ; ce que les rend capables de réaliser 


des performances très appropriées pour tout état de fonctionnement du réseau, [121]. 


2.2. Classification des dispositifs FACTS 


Depuis les premiers compensateurs, trois générations de dispositifs FACTS ont vu 
le jour. Elles se distinguent par la technologie des semi-conducteurs et des éléments de 
puissance utilisés [13, 14]. 

1. La première génération est basée sur les thyristors classiques. Ceux-ci sont 
généralement utilisés pour enclencher ou déclencher des composants afin de fournir 
ou absorber de la puissance réactive. Ils servent également au remplacement du 
changeur de prises en charge mécanique dans les transformateurs de réglage. 

2. La deuxième génération, dite avancée, est née avec l'avènement des semi- 
conducteurs de puissance réglables à la fermeture et à l'ouverture, comme le 
thyristor GTO. Ces éléments sont assemblés pour former des convertisseurs de 


tension ou de courant afin d'injecter des tensions contrôlables dans le réseau. 
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3. Une troisième génération de FACTS utilisant des composants hybrides et qui est 
adaptée à chaque cas. Contrairement aux deux premières générations, celle-ci 
n'utilise pas de dispositifs auxiliaires encombrants tels que des transformateurs pour 
le couplage avec le réseau. 

Les dispositifs FACTS peuvent être classés, selon leur mode de couplage ; 


parallèle, série et hybrides (série-parallèle), en trois catégories présentées ci-après : 


2.2.1. Dispositifs shunts. 


a. À thyristors, [12, 13, 14, 15]. 
e _ TCR (Thyristor Controlled Reactor). 
e __TSC (Thyristor Switched Capacitor). 
e SVC (Static Var Compensator). 
e _ TCBR (Thyristor Control Breaking Resistor). 
b. À thyristors GTO, [13, 15, 16]. 
e _ STATCOM (Static Synchronous Compensator)]. 


e  SSG (Static Synchronous Generator). 


2.2.2. Dispositifs série. 


a. À thyristors, [13]. 
e _ TCSR (Thyristor Controlled Series Reactor). 
e _ TCSC (Condensateur série commandé par thyristors). 
e __ TSSR (Thyristor Switched Series Reactor). 
e  TSSC (Thyristor Switched Series Capacitor). 
b. À thyristors GTO [13, 16]. 
e  SSSC (Static Synchronous Series Compensator). 


e _ GCSC(GTO Thyristor-Cotrolled Series Capacitors) 


2.2.3. Dispositifs hybrides. 
e  TCPAR ( Thyristor Controlled Phase Angle Regulator). 
e __ TCPST (Thyristor Controlled Phase Shifting Transformer). 
e _ TCVR (Thyristor Controlled Voltage Regulator). 
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e  IPEC (Interline Power Flow Controller) [15]. 
e  UPFC (Unified Power Flow Controller) [13]. 


2.3. Compensateur statique de puissance réactive (SVC). 


Le SVC est une installation shunt destinée à la production ou à l’absorption de 
l’énergie réactive, pour des objectifs de contrôle des paramètres spécifiques du réseau 
électrique, essentiellement la tension, [111]. 

D'un point de vue opérationnel, le SVC se comporte comme une réactance variable 
qui génère ou absorbe une puissance réactive afin de réguler la tension. La répartition 
adéquate des SVC sur le plan du réseau peut, dans une mesure importante, contribuer à 
l’amélioration de la capacité de transit et réduire les pertes de transmission [17]. 

Généralement, un SVC est constitué d'un ou plusieurs batteries de condensateurs 
fixes (FC) commutables soit par disjoncteur, ou bien par thyristors (Thyristor Switched 
Capacitor TSC) et d'un banc de réactances contrôlables par thyristors (Thyristor controlled 
Reactor TCR) et des filtres d'harmoniques. Le SVC est relié au réseau à travers un 
transformateur de couplage, [14], (Fig. 2.1.a). 

Les caractéristiques tension-courant d’un SVC sont montrées sur la figure 2.1, b. 
On y distingue trois zones de fonctionnement, [18] ; 

e L’intervalle V, < Vs < V2 définit une zone de réglage où l’énergie réactive 
est une résultante interactive entre le TCR et le TSC. La tension du nœud de 
connexion est, 

Vs = Vref + X 51.1 (2.1) 
La gamme de régulation du SVC est déterminée, alors, par les limites, 
-BCmax < B < Bimax 

e Pour V < V}(V,,) on définit une zone où le fonctionnement du SVC est 

purement capacitif, 


I 
Pr (2.2) 


e Pour V > V{(V,) la zone définie correspond au fonctionnement purement 
inductif du SVC est complètement inductif, le TCR est alors caractérisé par 


son énergie maximale ; les condensateurs état déconnectés. 
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Vs =- (2.3) 


Vs 


Shunt 
Transformer 


Ver 


Pente Xp 


+ 


Fig.2.1, a. Structure de base d'un SVC, b. 
Caractéristiques tension-courant, [19]. 


2.3.1. Réactance contrôlée par thyristors (TCR). 


Le TCR est composé d’une inductance en série avec deux thyristors tête-bêche (Fig. 
2.2, a). Les thyristors conduisent alternativement dans des intervalles dépendant de l'angle 
d'amorçage ” à ”. On définit l'angle d'allumage ” à ” à partir du passage par zéro dans le 
sens positif de la tension aux bornes du thyristor à allumer ; l'angle " © " étant l'angle de 
conduction des thyristors. 

Le système de contrôle les impulsions de gâchette des thyristors mesure soit la 
tension au point de connexion soit la puissance réactive dans la charge associée et calcule 
l’instant d’amorçage des valves. 

La pleine conduction des thyristors est obtenue pour un angle d'amorçage de 90°. Le 
courant est alors réactif et sinusoïdal, en retard de 90° par rapport à la tension du réseau. 
Lorsque a — 180°, on est en conduction nulle. La conduction partielle est obtenue pour un 
angle d'amorçage compris entre 90° et 180°, (Fig.2.2, b) [18, 20]. 

L’angle d'amorçage et l’angle de conduction s’expriment, dans leur relation, par, 

5=2{r-a) (2.4) 

L'augmentation de l'angle d'amorçage entraîne la diminution de l'amplitude du 

fondamental du courant. Cela se traduit par l'augmentation de la valeur de l'inductance du 


compensateur vue du réseau. 
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L'augmentation de l'angle d'amorçage rend la pollution harmonique du courant plus 


intense ; ce qui nécessite l’installation des filtres d’harmoniques. 


(a) (b) 


Thy2 


1 L ] 
D.T6 CRE CE] D#2 CET 


Fig.2. 2. Schéma d’une réactance contrôlée par thyristors. 


a. Schéma de principe, b. Tension et courant instantanés, [19]. 


2.3.2. Condensateur commuté par thyristors (TSC). 


Le TSC est composé d’un condensateur ” C " fixe branché en série avec une valve 
bidirectionnelle à thyristors (gradateur) et une bobine d’inductance d’atténuation " Le " 
(Fig.2.3, a). La commande envoie des impulsions sur les deux thyristors tant qu'il existe 
une demande de courant capacitif, donc de puissance réactive, (Fig.2.3, b). Le TSC 
fonctionne alors en pleine conduction durant cet intervalle de temps et il est bloqué si la 
demande en énergie réactive est nulle. Les thyristors sont utilisés comme des disjoncteurs 
électroniques et le condensateur n’est ainsi pas commandé en phase ; mais simplement 
enclenché et déclenché. L’inductance d’atténuation " L, " sert à limiter le courant en cas de 
fonctionnement anormal et à éviter la résonance avec le réseau. 

Pour avoir un minimum de perturbations transitoires, les instants de commutation 
sont choisis de façon à ce que la tension aux bornes des thyristors soit minimale. 
L’enclenchement est donc réalisé lorsque la tension résiduelle du condensateur est égale à 


la tension instantanée du réseau. 
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Fig.2.3. Schéma d'un Condensateur commuté par Thyristors. 


a. Schéma de principe, b. Courant et tension du condensateur [19]. 
2.4. Compensateur série commandé par thyristors (TCSC). 


Le TCSC est constitué d'une inductance en série avec un gradateur à thyristor, le 
tout en parallèle avec un condensateur (Fig.2.4, a). Le TCSC, placé en série dans une ligne 
de transport, permet de contrôler le flux de puissance et d'élever la capacité de transfert des 
lignes en agissant sur son impédance X7csc. 


X7cs 


A Région prohibée 


Cresc 


Thy Xrcsc (a) b) 


Fig.2.4, a. Structure de base d’un TCSC, b. Impédance du TCSC 


en fonction de l'angle d'amorçage (zone de fonctionnement, [19]). 


S1 les thyristors sont bloqués ” à =/80° ", le TCSC a une impédance fixe qui est 
celle du condensateur. Si les thyristors sont commandés en interrupteur électronique et en 
pleine conduction ” a = 90° ", l’impédance du TCSC est encore fixe et vaut l’impédance 
équivalente du condensateur en parallèle avec l’inductance. Si les thyristors sont 
partiellement ouverts " 90° < a <180° " (Fig.2.5), l’impédance du TCSC s’exprime, en 
fonction de l’angle d’amorçage "a", comme suit, 
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= Xc X1(0) 
Xrcsc (o)= AS 7 (2.5) 
ou bien, 
JjLw 
Xrcsc (a) = . (2.6) 
2 É — LS LCw? 
T 


a.- Thyristors bloqués. b.- Thyristors court-circuités c.- Conduction partielle. 


Fig.2.5. Régime de fonctionnement du TCSC. 


La figure 2.4, b illustre la variation du l’impédance du TCSC en fonction de l’angle 


" 


d’amorçage " a "des thyristors. On y distingue une zone où cette impédance n'est pas 
définie (zone de résonance). Dans le point de résonance, X:(@)= - Xc (Xrcsc(@)=) et la 
zone correspondante est prohibée pour le contrôle. En pratique, une varistance à oxyde 
métallique où MOV (Metal Oxide Vraistor) est insérée en parallèle afin de protéger le 


condensateur contre les surtensions [14]. 


2.5. Compensateur statique synchrone (STATCOM) 


Le STATCOM est défini comme étant un générateur synchrone fonctionnant, alors, 
comme un compensateur parallèle de l’énergie réactive ; dont le courant de sortie (capacitif 
ou inductif) peut être contrôlé indépendamment de la tension du réseau [15, 17]. Il est 
similaire au compensateur SVC sauf qu’il présente un temps de réponse plus court et il est 
sans inertie mécanique. 

Le STATCOM est principalement utilisé pour la compensation dynamique des 
réseaux, afin d’assurer la tenue de tension, d’accroître la stabilité en régime transitoire et 


d’amortir les oscillations de puissance. Son principe de fonctionnement est simple: un 
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onduleur de tension (convertisseur) produit, grâce à une source de tension continue, un jeu 
de tensions triphasées, en phase et couplées aux tensions d'une ligne via une inductance 
faible, correspondant en général à l'inductance de fuite du transformateur de couplage 
(Fig.2.6) [181]. 

L'échange d'énergie réactive avec le réseau se fait par le contrôle de la tension de 
sortie de l'onduleur Vsrar; laquelle est en phase avec la tension Vr du réseau. Le 
fonctionnement est similaire à celui du compensateur synchrone classique et peut être 


décrit de la façon suivante. 


La puissance réactive échangée avec le réseau est donnée par [21]: 


2 
V VeV 
OSsTAT = = R _ STAT cos(5 0) (2.7) 
STAT STAT 


avec 6, l’angle de la tension au jeu de barre sur lequel est couplée le STATCOM. 


Pour un STATCOM idéal (sans pertes actives), 8=0 , alors 


2 
VR 


= (Vr —Vsrar ) (2.8) 
STAT 


OSTAT = 


e Si Vsrar < VR, alors Qsrar devient positive, le courant circulant dans 
l’inductance est déphasé de + par rapport à la tension VR; ce qui donne un 


courant inductif, donc le STATCOM absorbe de la puissance réactive. (Fig.2.7, 
a). 
e Si Vsrar > VR, alors Osrar devient négative, le courant circulant dans 


l’inductance est déphasé de ee par rapport à la tension VR; ce qui donne un 


courant capacitif, donc le STATCOM fournie de la puissance réactive. (Fig.2.7, 
b). 
e Si Vsrar = VRr, le courant circulant dans l'inductance est nul et, par 


conséquent, il n’y a pas d’échange d’énergie. 
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Fig2..6. Schéma d’un STATCOM. a. Structure, 


b. Circuit équivalent, c. Echange de puissance. 


Is 
V V VRr Vs 
A 
a) courant inductif b) courant capacitif 


Fig. 2.7. Diagramme vectoriel de STATCOM. 


Le STATCOM présente plusieurs avantages ; 
e Un contrôle optimale de la tension, les phases sont contrôlées séparément 
pendant les perturbations du système ; 
e Bonne repense à faible tension : le STATCOM est capable de fournir son 
courant nominal, même lorsque la tension est presque nulle ; 


e Bonne repense dynamique c.à.d. le système repend instantanément. 


Un autre avantage du STATCOM réside dans le fais que l’énergie réactive 
échangée par ce compensateur peut être capacitive ou inductive à l’aide, uniquement, 
d’une inductance ; c.-à.-d, contrairement au SVC, il y à absence d’éléments capacitifs qui 
risque de provoquer des résonances avec des éléments inductifs du réseau [15]. La 


caractéristique statique de ce convertisseur est donnée par la figure 2.8. 
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Fig.2.8. Caractéristique du STATCOM. 


2.6. Compensateur statique série synchrone (SSSC) 


Le SSSC est un dispositif FACTS de la deuxième génération. Il est constitué d’un 
convertisseur de tension. Cette tension est induite à travers un transformateur dans son 
enroulement secondaire ; lequel est branché en série avec la ligne (Fig.2.9), autrement dit 
une source de tension connecté en série avec la ligne par l'intermédiaire d’un 
transformateur. Le SSSC agit, ainsi, sur le courant de la ligne en insérant une tension 
capacitive ou inductive de même fréquence. 

Si l’injection de la tension s’effectue avec un angle de 90° degré en arrière par 
rapport au courant de la ligne, on obtient une compensation série. Cette compensation est 
équivalente à un condensateur série à fréquence fondamentale. Si par contre, la tension est 
injectée avec un angle de 90° en avance par rapport au courant de la ligne, la compensation 
a le même effet que celui d’augmenter l’impédance réactive de la ligne. 

Le SSSC peut produire ou absorber la puissance réactive suivant la commande du 
convertisseur. Comparé avec le compensateur TCSC, le SSSC présente l’avantage de 
pouvoir maintenir la valeur de la tension insérée constante indépendamment du courant. Il 
est donc efficace tant pour des petites charges (faibles courants) que pour des grandes 
charges. Un autre avantage réside dans le fait que, contrairement au TCSC, le SSSC ne 
présente pas une zone morte à la résonance parallèle. De plus, sa taille est réduite en 
comparaison de celle d'un TCSC. Par contre, sa commande est plus complexe et son coût 


et plus élevé. 


34 


Mode capacitif 
" 


de 
Convertisseur contrôle 


Mode inductif 


Fig.2. 9. Compensateur synchrone statique série. 


a. Schéma, b. Caractéristique V-I 


2.7. Contrôleur de transit de puissance universel (UPFC). 


L’UPFC appelé aussi déphaseur régulateur universel (DRU) est l'équipement 
FACTS le plus puissant. Il est dit universel car 1l comporte toutes les stratégies de 
compensation de puissance. Il est capable de contrôler, simultanément et indépendamment, 
la tension, l'impédance, l'angle de transport de la ligne (déphasage) et d'intervenir dans le 
contrôle du transit de la puissance active. 

L’UPFC est constitué de deux onduleurs de tension. Un onduleur est connecté en 
parallèle au réseau par l'intermédiaire d'un transformateur shunt et l'autre est connecté en 
série avec le réseau. Les deux onduleurs sont interconnectés par une liaison à courant 
continu formé par un condensateur. Il s’agit en fait de la combinaison d’un STATCOM et 
un SSSC. Son schéma est représenté à la figure 2.10. L'onduleur 1 est utilisé pour 
satisfaire la demande en puissance active de l'onduleur 2 à travers la liaison à courant 
continu et 1l réalise aussi la fonction de compensation d'énergie réactive puisqu'il peut 
fournir ou absorber de la puissance réactive indépendamment de la puissance active qui le 
traverse. L’onduleur 2 injecte une tension contrôlable en amplitude et en phase et fournit 


les puissances active et réactive nécessaires à la compensation série [13,18]. 
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Fig.2.10. Schéma de principe de l'UPFC. 


2.8. Avantages et applications des FACTS. 


Les dispositifs FACTS, insérés dans un réseau électrique, peuvent toucher divers 
domaines tels que technique, économique et environnemental, et satisfaire plusieurs 
besoins. De manière générale, l'effet de ces équipements de compensation sur la 
performance des réseaux électriques se résume dans les points suivants ; 

e La contribution à l'exploitation optimale d’un système électrique par réduction 

des pertes de puissance et d’amélioration du profil de tension. 

e Meilleur contrôle et meilleure gestion de l'écoulement de puissance, de manière 

à réduire les surcharges dans les lignes ou les transformateurs. 

e Le maintien de la tension des nœuds à un niveau acceptable en fournissant (ou 

absorbant) de la puissance réactive appropriée. 

e Réduction de l’écoulement de la puissance réactive, ce qui permettra ainsi à la 

ligne de transporter plus de puissance active. 

e Possibilité de charger les lignes de transport près de leurs limites dynamiques. 

e _Amortissement des oscillations de puissances. 


e Amélioration de la qualité de l'énergie. 
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CHAPITRE 3 : REPARTITION OPTIMALE DES FLUX DE 
PUISSANCE EN PRESENCE DES FACTS. 


3.1. Introduction. 


Dans un sens général, l’optimisation des systèmes électriques, dans leurs états 
structurels et fonctionnels, est une problématique connue et, à juste titre très bien posée. 
Les résolutions des questions posées dans ce domaine sont diverses, selon les objectifs 
spécifiques fixés. 

Les moyens utilisés pour atteindre ces résolutions n’ont jamais cessé de se 
développer ; puisque les problèmes posés, eux-mêmes, sont sujets à évolution et à 
complexité incessantes. En effet, l’enjeu de l’énergie électrique, dont le besoin et la 
consommation ne cessent de s’affirmer, implique le développement pressant des systèmes 
électriques ; aussi bien sur le plan de la production que sur celui de la consommation. La 
complexité conséquente de leur structure et celle de leurs régimes dus à l’insertion de 
nouveaux types de charge (nouveaux centres de consommation) à caractéristiques, en règle 
générale, non linéaire, donc plus perturbatrices, justifie la recherche de nouveaux moyens 
et procédures, conformément adéquats, pour le traitement de cette problématique ainsi 
reposée et réactualisée. En règle générale, les procédures utilisées reposent sur le même 
formalisme analytique descriptif ; mais les moyens de résolution peuvent différer en 
fonction des apports progressifs en connaissances dans le domaine. 

Dans ce travail, nous nous basons sur ce formalisme descriptif, tiré de travaux qui 
ont précédé, essentiellement de [22+34], auxquels nous avons tenté d’appliquer une 


procédure propre à ce travail. 
3.2. Description analytique de l’écoulement de puissance. 


L’écoulement de puissance ou bien la répartition de charge (load flow or power 
flow) est d’impact très déterminant sur les paramètres du régime du réseau électrique. Son 
contrôle représente un moyen fondamental pour la planification et l’exploitation de ce 
dernier. Il permet, en fait, d’obtenir des informations sur le régime du réseau étudié et 
d’installer les conditions objectives d’une gestion adéquate dans diverses orientations telles 


que ; sécurité, fiabilité, économie, prévision. 
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En effet, à partir de données sur la consommation de l’énergie électrique et sur sa 
production spécifiées dans divers nœuds du système, la résolution du problème lié à 
l’écoulement de puissance permet d’atteindre les informations nécessaires sur l’état du 
processus électrique du réseau ; telles, les tensions nodales, les puissances de transits dans 
toutes les branches, les pertes dans les lignes, etc... [22]. 

Le traitement de l’écoulement de puissance consiste, fondamentalement et à travers 
une formulation mathématique adéquate, à déterminer les tensions complexes de tous les 
jeux de barres et, à partir desquelles, toutes les autres grandeurs ; telles que les courants de 
ligne, les puissances correspondantes, les pertes de puissances … 


Pour un réseau électrique à «n>» jeux de barres, le système d’équations des 


tensions nodales s’exprime par la relation matricielle, 


1 11 F, +. V 
rire @.D 
Le Y, Les 1. V, 
ou bien, tout simplement, par 
I bus — Vous *Vhus (3.2) 


où; 1h, = Vecteur de courant complexe injecté à chaque nœud, V,,, - Vecteur de la tension complexe de 
chaque nœud, Y,,, - Matrice d’admittance, elle contient les données concernant les lignes, les 


transformateurs et les condensateurs installés dans le système. 


Le courant /; s’écrit, 


A ME (33) 


avec, 


(3.4) 


La puissance injectée dans le système au jeu de barre i s’exprime comme, 


S;=V:*I; ; (3.5) 


soit, compte tenu de l’équation (3.3), 
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(0: —0 : -5.: 
paint ERP (3.6) 


Yi 


n 
S;=P-j0 = 
j=1 


où, 0; =0;-0;- Différence angulaire entre les jeux de barres à et j; 0;, 0 ;- Angles de la tension, 


respectivement, des jeux de barres i et j. 


Cette dernière équation montre que les puissances actives et réactives transmises 
dans un jeu de barre quelconque sont fonction du module et de l’argument des tensions aux 
niveaux des autres jeux de barres du système. En utilisant cette équation, on peut exprimer 


la puissance, respectivement, active et réactive, injectée dans un jeu de barre à sous la 


forme, 
n 
P.= Poi- Ppi = D VV;(e cos0, +b; sin0;) (3.7) 
j=1 
n 
Q; = QGi - On: = Y vv;(e; sin0;; —b;; cos0;) (3.8) 
j=1 
ou bien, 
n 
P= Po Pni = DV: Y; V; cos (b, 0-5; (3.9) 
j=1 
n 
Q; = Qi -Qni = > Vi lv; sin(6, —8; 5), (3.10) 
j=1 


où, gjÿet b; représentent, respectivement, la conductance et la susceptance de la ligne de transmission. 


L’écoulement des puissances dans une branche / du système peut être représenté, 
conformément au schéma de la figure 3.1, par les puissances apparentes d’extrémités 


comme suit, 


Li Si) ’ Sitli) j 
mt. Poe 
V V; 


Fig.3.1. Puissances de ligne. 


Sÿ=V;*Iÿ =] P?+0ÿ (3.11) 
Si =V;*I = P$ +0% (3.12) 


où, les puissances correspondantes actives et réactives ont pour forme, 


P; =V} 8; —V; y, ÿ ©0S 0; +bsin0;) (3.13) 
Pi =VPe Vi Vilejcos 0j +bjsin0;) (3.14) 
62 =_V?b; +V; ve, sin0;; —b;; cos 0) (3.15) 
Q; =_V?h; +V; vie; cos 0; —b ; cos 0;) (3.16) 


Les pertes correspondantes de puissance peuvent être, alors, déterminées comme 

suit ; 
SL, =Sÿ +S j (3.17) 
Dans les systèmes électriques, les nœuds sont caractérisés par une puissance 
complexe (S;) et un vecteur-tension (V;). Pour le traitement et calcul des régimes, il est 
nécessaire de disposer d’un minimum de données. Si on dispose de données sur la 
puissance uniquement, la tension, étant donc inconnue, sera à déterminer dans le processus 
de calcul. Si la tension est initialement donnée, la puissance inconnue sera alors à 


déterminer. A la base de ces grandeurs et de leur données relatives, on distingue trois de 


nœuds ; 


e _Nœuds P-V. Pour ce type de nœud, associé aux centrales de production, on spécifie 
la puissance active et le module de la tension. la phase de la tension et la puissance 
réactive sont inconnues et sont alors à déterminer après la résolution du problème. 

e Nœuds P-Q. Pour ce type de nœud, on associe généralement les charges. Ces 
dernières sont caractérisées par la consommation des puissances active et réactive. 
Les valeurs inconnues sont le module et la phase de la tension. 

e Nœuds V-6 : Pour ce type de nœud, on associe la centrale de production la plus 
puissante. Dans un nœud k (nœud de référence ou slack bus), on spécifie la phase et 
le module de la tension. Les valeurs à déterminer sont les puissances active et 


réactive. 


En raison de la nature non linéaire des équations de l'écoulement de puissance, 
plusieurs méthodes numériques ont été proposées pour résoudre ce problème. 
Généralement, les techniques les plus utilisées sont la méthode de Gauss-Seidel, Newton- 


Raphson et découplée rapide. 
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3.3. Résolution des équations de l'écoulement de puissance. 


La première méthode numérique utilisée était la méthode itérative de Gauss-Seidel. 
Cette méthode nécessite un grand nombre d’itérations pour les grands réseaux et un temps 
de convergence très grand. La méthode la plus connue est celle de Newton-Raphson ; 
laquelle représente un outil performant dans le traitement de la répartition des puissances. 
Cette dernière nécessite plus de temps par itération que celle de Gauss-Seidel, mais elle ne 
demande que quelques itérations pour atteindre la solution, même pour les grands réseaux, 
[23]. C’est la méthode (Newton-Raphson) adoptée et utilisée dans le cadre de ce travail de 
recherche, tout le long des divers traitements. 

Le principe de cette méthode est basé sur le développement en série de Taylor des 
équations des puissances (3.9) et (3.10). Comme les variations des puissances dans les 
deux équations sont liées aux variations de l’amplitude et de la phase des tensions, on peut 


écrire [24], 


 6P, GP, 

TR Er (3.18) 
j=17 le 7 
© 09; 20; 


En écrivant le système des équations précédentes sous la forme matricielle, on 


obtient, 
|A 64 © .. | 
OP... 0. OP. OP: | 
APy |_| 08: 08, OVy |xl A0» (3.20) 
AQ | |0Q . 9 0 .. 2] /|A" | 
08, 00, ©, OV, 
[A9] |29, ! 00, 0, 00, | lave 
| 06: 06, ô, ôV, | 


L’équation (3.20), peut être reformulée par l’équation suivante ; 


ol se G21) 
AQ!| |J3 Jal |AV 


, 


al 


avec ; J1,J3,J3,J4 - Respectivement, les sous-matrices formants la matrice Jacobienne qui 
représentent les dérivées partielles des P,Q par rapport à 6,V. Les éléments de la matrice 


Jacobienne évaluent, en fait, les dégres (sensibilités, vitesses) de variations des puissances 
actives et réactive par rapport au potentiel et à son déphasage. 


Les éléments diagonaux et non diagonaux de la sous matrice J, s’écrivent 


comme [22], 


Vi|sin(e,-0,-5;) (3.22) 


fl 


|A 


°F, _ 
m0, À 


n 
té 


Se 


œP, 


=}; 
00; 


V/sin(e,-e, -5;) (3.23) 


\A 


Analogiquement pour la sous-matrice J; ; 


©P; 2 

AR Yj|Vifeos be. = -55)+f;]hlcos (5) (3.24) 
_E 

©P, 

FA Yj Vilcos F5; 70; +0;) (3.25) 


Les éléments de J;, sont déterminés comme suit : 


ne > Val 1cos (:-0,-5;) (3.26) 
É 

40: 

BL ilbleos(e; -0, -5,) (3.27) 


De même, pour la sous matrice J4 ; 


09; L k 2 
D = Shin, -0,-8,)+1v;/frifsin(Gà) (3.28) 
A j=1 

i£j 

00: ; 

RL Visin(-5; -0,+0,) (4.29) 


Le calcul de l’écoulement de puissance par la méthode de Newton-Raphson se 


réalise suivant les étapes ; 


1. Préparation des données du réseau, qui incluent ; les données des lignes, des nœuds, 


des générateurs, des transformateurs, des charges, et l’erreur tolérable£, 
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2. Initialiser le compteur d’itérations, k = 0, 
Calcul de la matrice admittance Yys, 


4. Affecter les valeurs initiales aux grandeurs modules et phases des tensions ; 


habituellement, les valeurs initiales sont V: =1.0 p.u. et 89 =0 , 


5. Utilisez les estimations V° et 0° pour calculer la puissance active et réactive pour 
tous les nœuds P°, Q° en utilisant les équations (3.9) et (3.10), 
6. Calculer le vecteur des écarts de puissance active et réactive comme suit ; 


AP = psrec _ pb (3.30) 


AQ{° LE OPee _ o!° (3.3 1) 
Les écarts des puissances représentent la différence entre l’injection des puissances 


spécifiées et la valeur actuelle de l’injection des puissances calculées par (3.9) et 
(3.10). 

7. Utiliser les estimations V° et 6° pour calculer les éléments de la matrice 
Jacobéenne J,,J,,J3,J, en utilisant les équations de (3.22) à (3.29). 


8. Calcul de l’incrémentation des modules et des phases des tensions par 


l’inversement de la matrice Jacobéenne comme suit ; 


ab 2] ol 7 

9. Calculer les nouvelles estimations des modules et des phases des tensions selon les 
équations suivantes: 

VE VS HAVE (3.33) 


gk+1 2 pK + Agf (3.34) 


10. Si AP) <e et AO" <e l’algorithme converge, si non : 


11. Incrémenter k de et retourner à l’étape 5. 


3.4. Modèles des éléments du réseau. 
3.4.1. Générateurs 


En règle générale, les générateurs électriques, éléments du système d’énergie 
électrique et source de potentiel, sont modélisés par une force électromotrice débitant dans 


un nœud du système à travers une réactance. La valeur de ces grandeurs est fonction du 
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caractère du régime étudié (permanent, transitoire, subtransitoire) et du type de la machine 
(à pôles saillants, à pôles lisses). Les limites de leurs potentialités énergétiques sont 
déterminées par celles du courant inducteur et du courant induit (respectivement, du rotor 
et du stator) dans lesquelles limites s’installe la possibilité de contrôle de la tension et des 
puissances. Dans le cas considéré, le générateur est modélisé par une source de tension 
constante qui injecte, au niveau du nœud auquel il est connecté, une puissance active PG et 
réactive OG. 

La puissance générée par chaque unité est limitée dans un intervalle défini comme 
suit : 


LS PAS (3.35) 


ï 
06 < QG, SO (3.36) 

Les générateurs peuvent être donnés par deux modèles. Dans le premier, ils sont 
donnés par leur puissance Pi et leur potentiel de service V;; ils sont dits nœuds PV. La 
valeur de la puissance réactive est calculée de manière à maintenir la tension nodale à sa 
valeur de consigne. Si Qc; dépasse une des limites définies, alors elle sera fixée à cette 
même valeur et la tension ne peut plus être contrôlée. Le nœud se comportera comme étant 
un nœud PQ. 

Les générateurs peuvent êtres modélisés aussi par des puissances injectées 
constantes. Dans ce cas, les puissances 76; et Qc; sont données et la tension ne peut être 
contrôlée. Dès lors un générateur est connecté au nœud de référence (de la balance 
énergétique), tel que ses puissances active et réactive sont supposées infinies et sa phase 


sera choisie comme référence pour le système. La figure 3.2 illustre les deux modèles 


cités, [25]. 


Pi Qi 


Qc Qc 
Oimin 0 Oimax OGimin 0 QGimax 


Fig.3.2. Modèles du générateur. a- Symbole, b- Modèle 


PV classique, c- Puissance complexe constante. 
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3.4.2. Lignes. 


Une ligne est généralement représentée par son schéma électrique équivalent en 
«T » ; par son impédance série et par son admittance shunt ; soit, respectivement, 


Zi =ti+ Xi et Yi = gijo + Jbijo, (Fig.3.3). Généralement g << b et l’on peut considérer 
que g=0 et Yi = jbij Susceptance de la ligne, [25]. 


Zi Xi 


ri x 


ÿ 


ne 


Fig.3.3. Modèles en T de la ligne. 


Le nœud à de la ligne est en lien direct avec le nœud 0 (de potentiel nul) et le 
nœud j (de potentiel V;) par l’admittance, 


Y. = — + =Y..+ 3.37 
F :Z 2 + 2 ( ) 


1 
Fanny (3.38) 
] 2; J 


Analogiquement pour le nœud j, 


1 
F . =; + 
PR Do 2 


La matrice d’admittance nodale de liaison des nœuds et j est donc donnée par : 


Yÿ0 
Yÿ + > —}j 
Vous = > (3.39) 
Yi0 
di TT, 
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3.4.3. Transformateurs. 


Un transformateur simple à deux enroulements peut être représenté par le schéma 
de principe à liaison magnétique, (Fig.3.4, a) ou bien par son schéma électrique équivalent 
aux paramètres rapportés (ici primaire au secondaire) à travers le rapport de 
transformation ; soit a = Wy/W1, (Fig.3, b) et où, Lo, Zo — Paramètres, respectivement, 
courant et impédance de la branche d’aimantation, rapportés au secondaire. 


Pour ce schéma, on peut établir le système d’équations d’équilibre nodale ; 


où, les éléments de la matrice d’admittance peuvent être déterminés à partir de la mise en 


équation ou bien à partir de la mise en forme adéquate du schéma équivalent ; 


a. Mise en équation. 
Soit, donc, le système d’équations de l’équilibre du contour « 1-0-2 » et de celui du 


nœud «0 » ; 


aU, =al,Z, —-1,Z, +U, 
I 
1, =1,-— 
a 
En remplaçant le courant /: par sa valeur dans la première équation avec, à partir 
du schéma, 
Fe aU, -al,Z, 
0 A 
et, après quelques transformations, on obtient, 


Z 

f= : : DE “1 = 0 

éz(2)z * 
0 


De même en remplaçant 7 = a(1—1) dans la même équation avec, 
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I, = 
Pr CAE A 
ui 
W 
a. b. 


Fig.3.4. Représentation du transformateur. 
a.- liaison magnétique, b.- Liaison électrique. 


Ainsi, peut-on écrire ; 


Généralement l’impédance Z, est très grande, on peut supposer, donc, que Zo — 


æ, Ceci donne, 


où, en fait, &Z; + Z2 = Zr est l’impédance du transformateur, (Fig.3.5) ; 
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Fig.3.5. Schéma équivalent du 
TR avec liaison idéale 


ce qui donne ensuite, 


a al | 
In | Zr AU le a [Ui 
lue 
Zn Le a 
Ou bien, de manière générale, 
Me ue 
I; _ Y,; y; L/; = a a i 
lj Y; Y; U; Yi Y U; 
ä ÿ 
b. Transformation du schéma. 


La matrice des admittances 


peut être déterminée également à partir de la 
transformation du schéma de la forme « étoile » (Fig.3.4, b) à celle triangle (Fig.3.6). 


1 1;/a 


2,2 
aU;| 


Fig.3.6. Schéma transformé. 


Pour la forme triangle en peut définir, 


fe l 4 1 ÿ 1 
Zi Zn Zi 
Il 1 l 
ne Det lEl 
Zy Zn 
: 1 1 jé 
Fee + 
Zn Zoo Zu 


avec, 


2 
a Z.Z 
Zn =4 Z,+Z9 + 


2 


Z.,Z 
Ze Z 2 = À 


2; 
2 
VAE 
ARE PET NAT Es ee 
0 
Après substitution, on obtient, 
Z 
= 1 
2 2 Zo 
a" Z| —+1|+2Z, 
0 
1 
Y = 


Z; 
1 /a 1 Eu ; a, 
Ï Z Hal 
7 &Zz| 2+1|+2, 1  a=t+1] ? 
0 0 
ou bien, 
Z 
2e 
I, l a É "1 a A 
L Z 
I, feu) _q a | U, 
0 


Ainsi, pour le transformateur donné, la matrice d’admittances nodale s’exprime, 
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Ypus =| (3.40) 
3.4.4. Charges 


Les charges représentent les consommateurs de l’énergie électrique ; lesquels 
consommateurs peuvent aller de l’élémentaire à l’équivalent complexe très diversifié. En 
règle générale, une charge peut être modélisée par ; 

e Une impédance constante, 

e Un courant constant, 

e Une puissance complexe constante, 


e Une caractéristique statique complexe ; puissance complexe dépendante de la 


tension, S(U)= P(U)+ jO(U). 


Dans ce travail, la charge est donnée par une puissance complexe constante, 


indépendante de la tension du nœud auquel elle est branchée ; 


Si = Pi + JO (3.41) 
La puissance réactive O/,; peut être positive ou négative selon que la charge est de 


nature, respectivement, inductive ou capacitive, (Fig.3.7). 


Pa , d 


Fig. 3.7. Modèle de la charge. 


3.4.5. Eléments shunts. 


Pour des objectifs d’optimisation des régimes des systèmes électriques, on est tenu, 
très souvent, d’installer des éléments shunts pour divers contrôle, tel que des flux de 
circulation de la puissance réactive de ligne, de la tension, des pertes, de la capacité de 
transport... 


Pour le traitement des régimes, ces derniers sont introduits par leur admittance, 
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Yo = Gio + JBo. (3.42) 
Conformément à la nature dynamique des régimes du système, ces installations, en 
règle générale, sont de type réglable ; ils injectent au système de la puissance réactive aussi 


bien qu’ils en consomment, selon l’état du régime, (Fig.3.8). 


‘Système 


\ Electrique. 


Fig.3.8. Modèle de l'élément shunt ; 
Inductif- Capacitif. 


3.5. Modèles des dispositifs FACTS. 


Pour le contrôle des régimes du système aussi bien par les compensateurs classiques 
que par les FACTS, on peut recourir, en règle générale et selon l’objectif, à deux 
techniques ; l’une basée sur le contrôle par les paramètres du régime ; tels les puissances, 
le courant. Par cette procédure, on détermine les puissances réactives à injecter par les 


FACTS aux nœuds conformément aux critères d’optimisation donnés (Fig.3.9). 


S;=P;+ jO; S;=P;+jQ); 


Fig. 3.9. Modèle à puissance équivalente injectée. 
a. - Modèle avec FACTS ; b. - Modèle équivalent 
F.-Indice qui signifie FACTS. 
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Ces puissances équivalentes sont exprimées comme, 


SFACTS L pFACIS + jQFACTS (3.43) 


l 


TUE 0 (3.44) 


La matrice Jacobéenne est ainsi modifiée à travers les relations de sensibilité (degré 
de sensibilité) entre les paramètres du régime (AP, AQ, A, AU). La description de cette 
modification, suite à l’injection des puissances active et réactive par le dispositif FACTS, 


est donnée par la matrice Jacobéenne [26]; 


ML J' HE G45) 
AQ HP JAN AV 
avec : 
OP 
JV = Ji -——— 3.46 
I 1 ES (3.46) 
ap FACTS 
JE = T3 - —— 3.47 
2 2 ZV (3.47) 
FACTS 
Ps, (3.48) 
3 Oô 
FACTS 
PAR (3.49) 
4 OV 
L’écart de puissance active et réactive doit être modifié comme suit : 
AP = psPec Spie SP (3.50) 
: CIS 
AQ; = Qjÿ" - Of _0; (3.51) 
où, Pet QŸ*" sont, respectivement, la puissance active et réactive spécifiées au nœud i, 


PAS et QAR sont la puissance active et réactive injectées par un dispositif FACTS, 


P. et Q; sont respectivement la puissance active et réactive calculées par le programme 
d’écoulement de puissance. J,,J,,J3,J, sont les sous matrices de la matrice Jacobéenne 


sans dispositifs FACTS. 


La deuxième technique consiste en l’utilisation des paramètres du système ; tels les 


réactances. L’insertion des ces éléments dans la ligne change évidemment la structure de 


52 


dernière. Dans ce modèle, le dispositif FACTS vient, donc, modifier directement la matrice 
d'admittance de liaison nodale. Selon le type de FACTS modélisé, le dispositif peut être 


placé au milieu ou à une extrémité de la ligne (Fig.3.10), [27]. 


Fig.3.10. Insertion du dispositif FACTS dans une ligne. 


Les paramètres du modèle équivalent obtenu de la ligne sont déterminés et 


substitués à ceux de la ligne sans FACTS dans la matrice d'admittance nodale. Cette 


dernière est modifiée de la manière suivante, 


' ' ligne ligne FACIS FACIS 
_. Fe Yj | r Yi Y; 6.5) 
mod — | yr rs |T j : 
Y! y! ligne ligne FACTS FACTS ù 
J I : ji Y Ï Y Ji Y Ï 
avec : 
lign li 7ÿ0 
yligne  yligne Yi + Yÿ 
ll Ul _ 2 
Y ligne y ligne L V0 (3 5 3) 
Îi i =} Je : 


Selon le type de FACTS et sa position dans la ligne, seule une partie des 
coefficients de la matrice Y subit des modifications. Cette modélisation est celle utilisée 


dans le cadre de ce travail de recherche. 
3.5.1. Modèle du SVC à modification de la matrice d'admittance nodale. 


Le SVC est modélisé par une susceptance by variable, (Fig.3.11). Dans ce 
modèle, on considère la susceptance bc du SVC comme une variable de contrôle. La 


puissance réactive correspondante injectée dans le jeu de barres 1 est donnée par [28], 


2 
Osvc =—bsvc *V; (3.54) 
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Le signe « moins » indique que le SVC fournit de la puissance réactive au système 


lorsqu'il est capacitif alors qu’il en consomme lorsqu'il est inductif. 


line, = line 
Rj JXÿ 


J ’ijo 
Ysrc 


Fig. 3.11. Ligne avec SVC parallèle; variation discrète de la susceptance. 


Le courant du SVC en est déduit comme, 
Lsvc = j bsvc * V, (3.55) 
En présence du SVC dans le nœud ;, seul l'élément }; de la matrice d'admittance 
nodale est modifié ; l'admittance du SVC lui étant additionnée comme suit: 
Yi = Yi +Ysvc (3.56) 


Dans ce cas, la matrice d'admittance est modifiée de la manière suivante, 


Yÿ0 
. + PRE as FANNC . YE 
Yhus lys 1 |— . (3.57) 
Yi Yÿ Yÿ0 
MAP 
avec, 
Ysvc = J bsve (3.538) 


3.5.2. Modèle du TCSC à modification de la matrice d'admittance nodale. 
Le TCSC est un compensateur installé en série avec la ligne, utilisé principalement 


pour varier la réactance équivalente de la ligne, afin de contrôler le transit de puissance 


active et réactive dans le réseau. 
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- TCSC line |, - line 
EXy Rj"+jX; 


L/ 


Vio Vijo 


Fig. 3.12. Ligne avec SVC série; variation discrète de la susceptance. 


Le TCSC est donc représenté par une réactance variable X7Cc. La réactance 


résultante du modèle est, alors, la somme algébrique de la réactance du TCSC variable et la 
réactance réelle de la ligne ; 


line TCSC 
= A (3.59) 


La tension aux bornes du compensateur est définie par, 


Urcsc = j Xrcsc “1 (3.60) 
et la puissance réactive, par, 
Ufesc 
Qrcsc aus (3.61) 
TCSC 


La matrice d'admittance de la ligne est alors modifiée comme suit, 


! Yi0 ! 
a RE SU 
bus — e pale , jo (3.62) 
Ï “ji à Ï 
avec, 
! ! *D/ Il 
Yi = Gÿ + JB = (3.63) 


line + He ) 


F5 + JG 


3.6. Centrales Thermiques. 


Une centrale électrique est un centre de conversion énergétique destiné à la 
production de l’électricité. site industriel destiné à la production d'électricité. Les centrales 


électriques alimentent en électricité, à travers divers types de réseaux, divers 


consommateurs ; particuliers, industriels..., en règle générale, déployés et éloignés des 
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centrales. La production de l'électricité y est assurée par la conversion en énergie 
électrique d’une énergie primaire qui peut être de différente nature ; mécanique (force du 
vent, force de l’eau des rivières, des marées...), chimique (réactions d’oxydoréduction 
avec des combustibles, fossiles ou non, tels que la biomasse), nucléaire, solaire. 

Ces énergies primaires peuvent être renouvelables (biomasse), quasiment inépuisables 
(énergie solaire) ou des ressources dont la disponibilité est limitée dans le temps 
(combustibles fossiles). 

En fonction de ces énergies primaires, les centrales peuvent être de type thermiques 
classiques ou nucléaires, hydroélectriques ou à turbines à gaz, ainsi que des centrales 
décentralisées de l'énergie renouvelable solaire ou éolienne. 

La diversité des moyens de productions permet de faire face aux grandes 
fluctuations de la demande journalière, hebdomadaire et saisonnière d'une manière 
économique. Dans un contexte général, on distingue les centrales de base et les centrales 
de régularisation ou de contrôle. Les centrales thermiques assurent un rôle majeur dans 
l'équilibre d'un système électrique. En effet, les centrales thermiques nucléaires ou 
centrales thermiques à flamme (charbon, fioul, et gaz) assurent la demande de base [291]. 
Dans une centrale thermique classique, l’énergie électrique est produite, à partir de 
l’énergie thermique de combustion d’un fossile (charbon, pétrole, gaz naturel, etc.), 
préalablement convertie en énergie mécanique par une turbine à vapeur. Le schéma de 
principe d’un type simple de centrale thermique à vapeur est représenté par la figure 3.13. 

Le combustible, charbon, fioul ou gaz, est brûlé dans une chaudière. La quantité de 
chaleur dégagée transforme l’eau à vapeur. La vapeur produite est admise dans une turbine 
à vapeur où sa détente provoque la rotation des roues de la turbine, accouplée à un 


alternateur qui transforme l'énergie mécanique de la turbine en énergie électrique [29]. 


Chaudière Turbine Alternateur 
Entrée C T Sortie 
Energie (MW) 
| Electrique 
Fioul Eau 


Fig. 3.13. Centrale thermique à vapeur. 
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Les caractéristiques technico-économiques des centrales électriques sont 
déterminantes pour leur exploitation. Trois types de caractéristiques ont une influence pour 
l’exploitation d’une centrales électriques à court terme: son coût de production; ses 
contraintes techniques et sa fiabilité. Le plus important de ces trois caractéristiques est le 
coût de production [30]. Pour les centrales thermiques, le coût de production est lié 
principalement au coût du combustible défini par la modélisation des caractéristiques 


entrée-sortie des centrales. 
3.7. Modélisation des caractéristiques entrée-sortie des centrales thermiques. 


La caractéristique d’entrée-sortie d’une centrale thermique est définie comme la 
relation entre le taux de combustible en (Btu/h) ou son prix en ($/h) à l’entré et la 
puissance électrique en (MW) à la sortie de l’unité. Cette caractéristique est généralement 
représentée par un graphe de forme non linéaire, connu sous le nom de courbe de coût du 
combustible, (Fig.3.14). La forme de cette courbe est déterminée sur la base des données 
collectées à partir des essais effectués sur l’unité génératrice. La fonction associée à cette 


courbe est approximativement présentée sous la forme quadratique, [31]: 
F(Pa)= 4; +b, Pa +c Fe ($/h) (3.64) 


où, P,, est la quantité produite ; a,,b, ,c; sont les pondérations du coût de l’unité 


génératrice. 


P;"i Puissance P;"ax 


Fig. 3.14. Courbe modèle du coût du combustible. 


57 


L’appel en puissance utile active par le consommateur se traduit par l’augmentation 
du couple de rotation de la turbine ; laquelle augmentation implique la nécessité d’un débit 
plus grand du combustible. Tout ce processus s’accompagne, dans un général, par une 
variation du coût de la production de l’énergie électrique. Cette relation entre coût de 
production de l’énergie et son débit utile est d’autant plus complexe que les centrales à 


traiter sont plus diverses (nombre, type, caractéristiques,..….). 
3.8. Problème de l’écoulement optimal de puissance. 

3.8.1. Formulation du problème. 

L’écoulement optimal de puissance permet d’optimiser un ensemble de variables de 
contrôle afin de minimiser une fonction objectif prédéfinie tout en satisfaisant des 


contraintes d'égalité et d'inégalité. Mathématiquement, OPF peux être formulé comme un 


problème d’optimisation non linéaire sous contraintes de la forme [36] ; 


Minimiser F(xu) (3.65) 
Sujetà: g(xu)=0 (3.66) 
Et A(xu)<0 3.67) 


avec ; F(xu) - Fonction objectif, gx u) - Représente les contraintes d'égalités, hxu) - 
Représente les contraintes d'inégalités, u - Vecteur des variables de contrôle, x - Vecteur 


des variables d'état. 
3.8.2. Variables. 


Variables de contrôle. Les variables de contrôle, appelées également variables de 
décision, représentent l'ensemble des paramètres du système dont la valeur peut être 
ajustée ou contrôlé dans le processus d'optimisation pour obtenir les solutions optimales ; 
ces variables peuvent inclure : 


" Pc : Les puissances actives générées par les alternateurs sauf celle de référence. 


“  VG : Les modules des tensions des nœuds de contrôle (PV) et de référence. 
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“" T :Les rapports de transformation. 
"Qc : Les puissances réactives des compensateurs shunts. 
"  Xxcsc : La réactance du TCSC. 
" _ Bsyc : La susceptance du SVC. 
Par conséquent, le vecteur des variables de contrôle du système « u » est défini par, 


T 
= Po, Peyer Ver Vence» Qi Qcue Ti TT 


(3.68) 
X rcsc j: à X rcsc NTCSC ? Bsvc, Bsvc NSVC 


où, NG, NC, NT, NTCSC , NSVC sont, respectivement, le nombre de générateurs, le nombre 
de compensateurs shunts, le nombre de transformateurs, le nombre de dispositifs TCSC et 


le nombre de dispositifs SVC. 


Variables d’état. Ce sont les variables qui décrivent les changements d’état du 
système impliqués par les variables de contrôle. Elles peuvent inclure, 

" Pc, : La puissance active générée du jeu de barres de référence. 

" _V; : Les modules des tensions des jeux de barres de charge. 

" QG : Les puissances réactives des générateurs. 

" _S, : Les puissances apparentes transitant par les lignes de transmission. 
Par conséquent, le vecteur des variables d’état du système « x » est défini par : 

“A =|?,, Va Vans QG" QGye SL "Si | (3.69) 
où ; LN, nl sont, respectivement, le nombre des jeux de barres de charges (PQ) et le 
nombre de lignes de transmission. 

Ces variables ne sont pas directement contrôlées dans le processus d'optimisation. Ce sont 


des variables inconnues et généralement sont obtenues par la résolution de l'équation de 


l'écoulement de puissance [29]. 
3.8.3. Contraintes. 


Contraintes d’égalités. Ces contraintes sont l'image des lois physiques gouvernant 


le système électrique ; elles représentent les équations d’écoulement de puissance, dont ils 


59 


expriment l’égalité entre la production et la consommation de la puissance active et 


réactive à chaque nœud du réseau. Ces contraintes sont données par, 


j=1 
NB 
j=1 


Contraintes d’inégalités. Ces contraintes représentent les limitations physiques et 
techniques des différents équipements installés au réseau électrique ; elles sont résumées 
comme suit : 

e Contraintes du générateur. 
La tension du générateur, et les puissances actives et réactives fournies sont limitées 


par leurs limites inférieures et supérieures comme, 


VE Ve SE dela NG (3.72) 
PSP = LR NC (3.73) 
OP < Qc, <QP, i=1,..,NG (3.74) 


e Contraintes du transformateur. 


Le rapport du régleur en charge est limité par un minimum et un maximum comme, 


HET SA NT (3.75) 
e Contraintes du compensateur shunt. 


Les puissances réactives des compensateurs shunts doivent se situer dans les limites, 


QU < Qc, <QP, i=1...,NG (3.76) 


e Contraintes du compensateur statique d’énergie réactive (SVC). 


La susceptance du SVC est limitée comme suit : 


be, <bsvc, bye 1= 1. NSVC (3.77) 


e Contraintes du condensateur série commandé par thyristor (TCSC). 


La réactance du TCSC doit répondre, également, aux limites, 
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x Ken. Hi NTGSC (3.78) 


TC, - TCSC;? 

e Contraintes de sécurité. 
Le réseau est dit en état de sécurité s’il respecte les contraintes de sécurité suivantes, 

Ÿ”_ Puissance des lignes. 

Pour maintenir le système électrique en état de sécurité, les lignes de transport et les 
transformateurs de puissances doivent répondre à des limites des puissances de transit. Ces 
limites ont une relation directe avec les effets thermiques et électrodynamiques engendrés 
par le transit des puissances (passage du courant) et elles sont strictement déterminées par 
les « tenues thermique et électrodynamique » des équipements de force des réseaux. Dans 
ce qui suit, elles sont introduites par une contrainte d'inégalité, qui limitera l’écoulement de 


puissance sur chaque branche par une valeur maximal comme suit, 


Si, Dar i=1,..,NL (3.79) 
où ; $, - Puissance apparente circulant dans la branche / , définie par les équations (3.15) et 
(3.16), S “hi - Puissance maximale admissible dans la branche /, correspondante à la valeur 


maximale 1°" du courant de transit, déterminé par une tenue thermique admissible de 
l 


l’équipement de force. 


YŸ” Tension des nœuds. 
Pour des raisons de qualité, de sécurité et de stabilité du système électrique, les tensions ne 


doivent pas dépassées leurs limites admissibles données par, 


VV eV ent (3.80) 


L'écart de ces limites peut prendre des valeurs diverses selon la position hiérarchique du 
réseau et les exigences en qualité de l’énergie électrique (pratiquement jusqu’à 10% de la 


valeur nominale). 


3.8.4. Manipulation des contraintes. 
Il se peut, durant l’évaluation du problème d’optimisation, que certaines solutions 


violent une ou plusieurs contraintes de sécurité ; ces dernières ne peuvent pas être 
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considérées comme une solution optimale au problème d’optimisation, même si elle 
optimise la fonction objective. Il existe plusieurs techniques pour remédier à ce problème ; 
parmi lesquelles on a choisie celle qui pénalise les solutions violant les contraintes de 
sécurité. Le principe de base de cette méthode consiste à modifier la fonction objective en 
lui ajoutant une fonction de pénalisation ; la nouvelle fonction objective aura alors pour 


forme [22], 


Fcib (xu )= F(xu )+ Fpénatté (3.81) 


où, Fc (2 ) est la fonction globale proposée, Fpgnaié est la fonction de pénalité. 

Les variables de contrôle ne sont pas considérées dans la fonction objective globale, 
elles sont maintenues automatiquement à leurs limites par le programme d’optimisation. 

Les contraintes d'inégalité des variables d’état, qui contiennent la puissance active 
générée du jeu de barres de référence, les tensions du nœuds des charges, les puissances 
réactives des générateurs et les puissances apparentes transitant par les lignes de 
transmission sont multipliées par des facteurs de pénalité et inclues dans la fonction 
objective globale sous forme quadratique afin de maintenir leurs valeurs proches de leurs 


limites de fonctionnement. En incorporant toutes ces contraintes, on aboutit à la formule 


suivante, 
À NL “is NG im 2 
Fotb = F(xu)+Ap Cr pi) +ky À (vi, vi) + D (Oo 0) 
: Êi Fe (3.82) 
Moose) 
i=0 
avec, 


Rp (re, —- pain) : Fonction de pénalisation de puissance active du nœud de référence 


NG 
19 S (0 - oi) : Fonction de pénalisation des puissances réactives des générateurs 
i=1 


nl 
ks S(s, —- sim ) : Fonction de pénalisation des puissances transitant par les lignes. 
i=0 


À P : Facteur de pénalisation de la violation de la puissance active du nœud de référence. 
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À : Facteur de pénalisation de la violation des tensions des nœuds de charge. 

ÀQ : Facteur de pénalisation de la violation des puissances réactives des générateurs. 

1 : Facteur de pénalisation de la violation des puissances transitant par les lignes de 
transmission. 


Dans la fonction objective globale les termes 21m, vin Qlm slim sont définis 
J 8 G L; G; li 


dans les équations suivantes : 


min . min 
F6, ST F6 < le, 
lim _ max : max 
F6 = F6, si FG > F6, (3.83) 
. min max 
en ST F6, < EG <F6, 
A " V.< por 
pin y AS (3.84) 
V _ par <V,< pre 
0" Fe Qc, < Op" 
lim _ max . max 
DM sm 05>0P (3.85) 
. min max 
QG; s 6, < QG; < O6, 
Si si Si, eo 
L L 
SHSISER Ge Sn (3.86) 
Si. Si Serres < SX 
i l; i l; 


3.8.5. Fonctions objectives. 


Le problème de l’écoulement optimal de puissance peut être classé, selon le type de 
la fonction objectif, en deux sous problèmes ; l’écoulement optimal de puissance active 
(Optimal active power flow OPF), également appelée dispatching économique (economic 
dispatch DE) et l’écoulement optimal de puissance réactive (Optimal réactive power flow 
ORPEF) (Fig.3.15). Ces deux problèmes constituent l’objet d’intérêt dans l’étude qui suit ; 


orientée vers l’optimisation des fonctions objectives suivantes [32], 
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Contraintes d'épée OPF — Contraintes d’inégalités 


Ecoulement Optimal Ecoulement Optimal 
de Puissance Active de Puissance réactive 


Fonctions Objectif 


sécurité du 
système 


stabilité 
de tension 


coût de 


pertes actives 
et réactive 


Multi Objective OPF 


Fig.3.15. Basic structure of multi objective OPF strategy. 


déviation de 
tension 


combustible 


a. Répartition Optimale des Puissances Actives. 
Fonction objective 1: Minimisation du coût de combustible. 


La minimisation du coût du combustible est l'objectif le plus posé dans la 
formulation du problème d'OPF ; le but étant d’exploiter le réseau électrique avec un 
minimum de coût. Le problème consiste à trouver les puissances actives optimales à 
générer par les centrales de production auxquelles correspondent les frais minimaux de 


production. Cette proposition se traduit par la minimisation de la fonction suivante [33], 


NG NG 
F=Ÿ f; = Y'a;+b; PGi+ci Pi (S/h) (3.87) 
il 


i=1l 


où, f, est la fonction du coût du combustible pour chaque centrale électrique. 


b. Répartition Optimale des Puissances réactives. 
L’écoulement optimal de la puissance réactive s’exprime en la répartition optimale 


des flux de cette dernière dans la structure des réseaux électriques. L’objectif est de 
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minimiser les pertes de transmission, le transit de puissance, tout en garantissant la sécurité 
et la stabilité de système. Cette proposition se traduit par la minimisation des fonctions 
objectives suivantes. 

Fonction objective 2: Amélioration du profil de tension. 

La tension du jeu de barres est l’un des indices de qualité et de sécurité les plus 
importants dans les systèmes d'énergie. L’amélioration du profil de tension est obtenue si 
la tension du jeu de barres de charge est maintenue dans un écart admissible au voisinage 
de sa valeur spécifiée, généralement, nominale (1.0 pu). Dans ce cas, la fonction objective 
est minimisée conformément à cette proposition, où, la déviation DV de la tension des jeux 


de barres de charge, est définie par [29]: 


NL NL 
pv=>|v-vr# =Y|v-10| (3.88) 
i=1 i=1 


où, v' - Module de la tension de référence. 


Fonction objective 3: Amélioration de la stabilité de tension. 

La stabilité de tension est définie comme étant la capacité du réseau électrique à 
maintenir la tension des nœuds dans les limites permises, dans les conditions de 
fonctionnement normales ou suite à une perturbation quelconque [34]. D'autre part, Un 
système entre dans un état d'instabilité lorsqu'une perturbation provoque une décroissance 
progressive et incontrôlable du plan du potentiel du système. On assiste, alors à un 
effondrement total du potentiel. 

L'impact du plan du potentiel sur les caractéristiques de fonctionnement des 
systèmes électriques est très déterminant. De ce fait, de nombreuses recherches ont été 
orientées pour le traitement et la maîtrise des processus de contrôle de la stabilité du 
potentiel. On en cite les travaux de Kessel et Glavitsch [35] ; dans lesquels ces derniers 
ont développé des analyses et proposé un indice de stabilité de la tension dit « L-index ». 
Cet indice donne une image complète de la tension du nœud. Les valeurs données varient 
entre 0 et 1 ; la valeur 0 indique qu’aucune charge n’est connectée, tandis que la valeur 1 
indique la présence d’un écroulement de la tension. En effet, le jeu de barre avec la valeur 
de « L-index » la plus haute sera le jeu de barre le plus critique dans le système. Donc 
l’amélioration de la stabilité de la tension est faite par la minimisation de l'indice de 
stabilité. 

L-index est déterminé à partir de l'équation de base de l'écoulement de puissance 


déjà discuté précédemment, et qui est formulé comme suit : 
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Thus = Vous Vs (3.89) 
En séparant les jeux de barres de charge (nœuds PQ) des jeux de barres de 


génération (nœuds PV), la formule précédente peut être réécrite de la façon suivante: 


Da ARR 71117 
ME FE h (3.90) 


OÙ, Yu, Ma-YorYca SOnt des sous matrice de}, ; V,,1, Sont respectivement les tensions 


et les courants des jeux de barres de charge ; V,,14 Sont respectivement les tensions et les 


courants des jeux de barres de génération. 


On peut écrire la relation précédente avec une autre formulation comme suit, [36]: 

VE IL | | Hi Hi || 11 

1G Vo] LHcz Hcc]lVG 
où, H est déterminées par l'inversion partielle de la matrice admittance}},, ; 
HyHic-Hcr.Hca SOnt des sous matrice de H. 


La matrice H est donnée par : 


pe] 71 216 | (3.92) 
YGozZiz YGc -Yc1Zi ic 


avec, 
Zi =Y (3.93) 
Finalement, l'indice L des jeux de barres de charge j noté par L; est formulé comme 
suit: 
NG v. 
= TA ,  j=L2....NL (3.94) 


Pour assurer la stabilité en tension la condition L; <1 ne doit pas être violée pour 


tous les jeux de barres du système. 


L'indice de stabilité de tension L-index (Z 


ne) 


peut être donné par l’équation, 


L- index = Lax =max(L;),  j=12...NL (3.95) 


î 
Donc l’amélioration de la stabilité de la tension est faite par la minimisation de 

l'indice de stabilité L pour chaque jeu de barres de charge et ensuite la minimisation de 

l'indice maximal global L,.. Donc la fonction objective est donnée par l’expression, 


min L—index = min (max (z, ) (3.96) 
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Fonction objective 4: minimisation des pertes actives de transmission. 
Dans ce cas, la fonction objective à minimiser représente les pertes actives totales 


dissipées dans les lignes et qui peut formulée par l’expression suivante, 


NB NB NB NB 
F=YP = Pa Pi => V;(G; COS (,-0,)-8; sin (6, -0;)| (3.97) 
i=1 i=1 i=1 j=1 


Fonction objective 5: minimisation des pertes réactives de transmission. 


La fonction objective représente des pertes réactives totales de lignes ; elle est 


donnée par l’expression, 


NB NB NB NB 
F=-Y 0) O0) On =V:D.V;fGysin (6; -0;)-B; cos (6; -6;) (3.98) 
i=l i=l i=l j=1 


Fonction objective 6: maximisation de la sécurité du système. 

La fonction objective, dans ce cas, consiste à contrôler les transits de puissances de 
manière à réduire, voire supprimé, les surcharges dans les lignes et dans les 
transformateurs. L’expression à minimiser les transits de puissances dans les lignes est de 


la forme, 


Û (3.99) 


avec ; $;°"" - Puissance maximale admissible dans la branche /, $; : Puissance apparente 
l 


circulant dans la branche 1, donnée par les équations (3.15) et 3.16), P,,0;- 


Respectivement puissance active et réactive circulant dans la branche /, données par les 


équations (3.11-14). 
3.9. Résumé de la formulation de l'OPF (problème complet). 


NL 
ne 
Four) = Fi33456 000 +p (Pa Pa] ose 


:. |2 
lim 
. v, va) 
i=1 


L 
(3.100) 


nl 


NG as D 

Am m 

sis D loanoin uns disent) 
i=1 


i=0 
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NB 

1 

/ NB 
j=1 


Var eV SV 
PSS 
0" < Qc; < Our 
Ta < T, < mo 
2" < Qcc, < an 
bye, <bsve, <bSve. 


min max 
< < 
XrcsC; £ Xrcsc; < X7csC, 


S,, < sax 
avec : 

Si = JP2+0? 

P; = V8; vi; ( ÿ COS 0; + b;; sin 0ÿ) 


2 . 
Pa = Vis VV ie j cos ; +bj sin0 ;) (4.101) 
Qj = Vb; +VV (es sinO;; — b;; COS 0) 


2 
O ; = —V,; b;; +VV(e ; COS 6;; —b ; COS 0) 
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CHAPITRE 4 : METHODES D'OPTIMISATION ET PROPOSITION. 


4.1. Problèmatique et définitions. 


Un problème d’optimisation est défini, au sens général, par un ensemble de 
variables « x », une fonction objective « f » et un ensemble de contraintes d'égalité et/ou 
d'inégalité «C ». L’ensemble des solutions possibles du problème forme l’espace de 
recherche «S ». Ce dernier contient un sous-ensemble x € $S qui représente des solutions 
réalisables respectant les contraintes «C » du problème ; à chaque solution «x » est 
associée une valeur «f(x)» qui représente sa qualité. La résolution du problème 
d'optimisation consiste à explorer l’espace de recherche «S » afin de trouver une solution 
réalisable x* € S qui minimise ou maximise la valeur de « f(x) » [50]. 

Quelque soit le type du problème d’optimisation, ce dernier est défini par le 6-uplet 
<D, f,S, C, x, mode>. Où « D » représente les données du problème, « f » est la fonction à 
optimiser (à maximiser au minimiser), appelée aussi fonction coût ou fonction objectif, qui 
associe à chaque solution «x » une valeur numérique « f(x) » représentant la qualité de 
«x »; «S » représente l’ensemble des solutions possibles du problème traité, usuellement 
appelé espace de recherche ou espace des variables. «C » représente les contraintes 
d'égalité et/ou d'inégalité que doit satisfaire une solution afin d’être admissible. « x » est un 
vecteur de variables dites de décision qui appartient à l’espace de recherche &«S » et qui 
représente les solutions réalisables (admissibles). « Mode » indique le type du problème ; il 
permet de savoir est-ce qu’on doit minimiser ou maximiser les valeurs des solutions de 
«x » [52]. 

Un problème d’optimisation est usuellement formulè comme un problème de 


minimisation et est écrit sous la forme : 


min f(x) 

telque 
gi(x)< O,i=1.., m, (4.1) 
hG)=0 215 

XEScCR 


où, «x» est le vecteur des variables de décisions. «S » représente l’ensemble des 


solutions réalisables et f(x est la fonction objective ; « h;,g; » représentent les contraintes, 
respectivement, d'égalités et d'inégalités. 
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Nous présentons, dans ce qui suit, quelques définitions de base liés au problème 


d'optimisation [52, 53]. 


e Les contraintes : sont des restrictions imposées par la nature et les caractéristiques 
du problème sur les solutions proposées. 

e  L’espace de recherche : est composé de l’ensemble de valeurs pouvant être prises 
par les variables. 

e Espace des objectifs : ensemble image de l’espace de recherche occupé par les 
fonctions objectives. 

e solutions réalisable : un point X, est appelè un point admissible si : 
XA€S et si les contraintes d'égalité #;(x)=0 et d'inégalité g;(x)<0 sont satisfaites. 

e Domaine non-réalisable : région de l’espace où les contraintes sont violées. 

e Le voisinage : Le voisinage V(x) d’une solution «x » est un sous ensemble de 
«S », dont les membres sont des solutions proches (voisines) de la solution « x ». 
En effet, une solution « x* » est considérée comme une voisine de « x », si elle est 


le résultat d’une modification légère autour de la solution « x ». 
° L’optimum local: X de l’espace de recherche «S » est un minimum local de la 
fonction « f », s’il existe un voisinage de « X * » noté « v(x*) », tel que: 
vXev(x"):f(xX")<f(X) (4.2) 
e  L’optimum global: X°eS est un minimum global de la fonction « f » si: 
vXes:f(X )<f(x) (4.3) 


e Une fonction multimodale présente plusieurs minima (locaux et globaux), et une 


fonction uni modale n’a qu’un minimum, le minimum global. 


La figure 4.1 montre une fonction multimodale à deux optima locaux et un optimum 


global [54]. 
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Fig.4.1. Fonction multimodale à deux optima locaux et un optimum global. 


Un problème d’optimisation peut être statique ou dynamique (dans lequel la 
fonction objectif change avec le temps), mono ou multi-objectifs (à plusieurs fonctions 
objectifs souvent contradictoires) et avec ou sans contraintes. En outre, un problème 
d'optimisation peut être à variables continues (les variables composant une solution 
donnée sont de type réel), discrètes (variables de décision de type entier), ou combinatoire 
(mixte discret et continu). Un problème d’optimisation peut être aussi de type linéaire si la 
fonction objective est linéaire par rapport à toutes les variables et non linéaire dans le cas 
contraire. Ainsi, un problème d’optimisation peut être, par exemple, à la fois combinatoire 


et multi-objectifs [55]. 
4.2. Optimisation multi-objective : résolutions et stratégies. 

Un problème d’optimisation est dit multi-objectif ou multicritère lorsqu’1l nécessite 
l’optimisation simultanée de plusieurs fonctions décrites par les mêmes variables de 
décision. La forme générale d’un problème d’optimisation multi-objectif est donnée par le 


système d’équations suivant : 


Optimiser f(9 = (f(x, 20, fn(9)x e S (4.4) 
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avec n22 est le nombre de fonctions objectifs que l’on cherche à optimiser ; dans le cas où 
n=1, le problème d’optimisation multi-objectifs est réduit à un problème d’optimisation 
mono-objectif. 

Les problèmes multi-objectifs ont la particularité d'être beaucoup plus difficiles à 
traiter que leur équivalent mono-objectif. La difficulté réside dans la présence de conflits 
entre les diverses fonctions objectifs ; puisque les solutions optimales des unes ne 
correspondent pas généralement à celles des autres fonctions, cela signifie que 
l’amélioration d’une fonction entraîne la détérioration de l’autre ou des autres. En effet, on 
n’est plus ici à la recherche d’un unique optimum global (le cas d’un problème mono- 
objectif), mais plutôt d’un ensemble de solutions qui offrent un bon compromis entre les 
différents objectifs. Les méthodes de résolution des problèmes multi-objectifs sont donc 
des méthodes d’aide à la décision car le choix final sera laissé au décideur [55]. 

Pour la résolution d’un problème d’optimisation multi-objective, la communauté 
scientifique a adopté deux types de comportement qui traitent les objectifs simultanément 


ou séparément comme suit [58] : 


e Transformer le problème multi-objectif à un problème mono-objectif, en utilisant 
généralement des méthodes d'agrégation. 
e Utiliser directement la notion de dominance au sens de Pareto dans la sélection des 


solutions générées. 


4.2.1. Méthodes de transformation : multi-objectif en mono-objectif. 


Cette classe de méthodes, transformant le problème multi-objectif en un problème 
mono-objectif, génère donc une seule solution et non pas un ensemble de solutions Pareto 
optimales. La conversion est habituellement faite en agrégeant tous les objectifs dans une 
fonction pondérée. Parmi ces méthodes, on trouve principalement les méthodes 
d’agrégations (additives et multiplicatives), la méthode de pondération, l’approche min- 
max, la méthode de contrainte, l’approche du vecteur cible, et la méthode de 
programmation par but... C’est la méthode de pondération adoptée et utilisée dans le cadre 


de ce travail de recherche [591]. 


4.2.2. Méthode de pondération, (weighting method) 
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L’idée de cette méthode est d'associer à chaque fonction objectif un coefficient de 
pondération et à faire la somme des fonctions objectifs pondérées pour obtenir une 
nouvelle et unique fonction objectif. Cela modifie le problème multi-objectif, comme 


défini dans l’équation (4.4), en un problème simple objectif de la forme [58]: 


Optimiser f(x) = w.f1(0) +w2.f2(0)+uw fn, x ES (4.5) 


où », représentent des coefficients de pondération appropriée. w, doivent être choisis en 


fonction des préférences associées aux objectifs, ce qui est une tâche difficile. Cette 


méthode est simple à mettre en œuvre et elle est d’une grande efficacité. 


4.3. Notions de dominance et d’optimalité de Pareto. 


L'idée d’utiliser la dominance au sens de Pareto pour résoudre les problèmes 
d'optimisation multi-objective, a été proposée par le mathématicien Italien Vilfredo Pareto 
[34]. Le principal avantage de cette approche, c’est l’optimisation simultanée d’objectifs 
contradictoires. En effet, Une solution appartient à l’ensemble de Pareto si et seulement si, 
il n’existe pas une autre solution qui améliore au moins une fonction et ne dégrade pas les 


autres fonctions objectifs. Mathématiquement le problème peut être définit comme suit : 


Définition : Une solution réalisable xeS est Pareto optimale (dominant ou non dominée) 
si et seulement si il n’existe pas de solution x" es telle que x domine x. 
On dit que le vecteur y=[1,»2.., y] domine le vecteur z=[2,2..., z], dans le cas d’une 


minimisation d’objectifs, si et seulement si [57]: 


vie{L2,..k} fiO< fn ef2,.. k}: F0) < F0) (4.6) 
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Fig. 4.2. Front Pareto pour un problème d'optimisation à deux fonction objectif. 


Ensemble de solutions Pareto optimales et frontière. Une solution «x » est 
considérée comme optimale au sens de Pareto, ou Pareto-optimale s'il n'existe aucun point 
de l'espace faisable « S » qui la domine. Ces points sont également appelés l'ensemble de 


Pareto ou également l'ensemble des compromis optimaux [56]. L'ensemble de solutions 


Pareto PO* optimales est donné par, 
PO hesta ES: fR)< FU) (4.7) 


L'ensemble de tous les points Pareto-optimaux forment la frontière Pareto- 


optimale. Le Front Pareto PF* est défini comme suit (Fig.4.2) : 
PF =k=F=(f(00 fre P} (4.8) 


4.4. Méthodes d'optimisation. 


Une grande variété de techniques d'optimisation a été appliquée pour résoudre des 
problèmes d’optimisations. Ces méthodes peuvent être classifiées en deux groupes : 
déterministes et stochastiques. Une méthode d'optimisation est dite déterministe lorsqu’elle 
conduit, pour une solution initiale donnée, au même résultat final. Ce sont en général des 
méthodes exactes, efficaces, mais qui nécessitent une configuration initiale (point de 


départ) pour résoudre le problème. Les méthodes stochastiques sont la classe générale des 
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algorithmes et des techniques qui emploient des mécanismes de transition aléatoire pour 
trouver des solutions optimales à des problèmes difficiles. Cette caractéristique indique que 


plusieurs exécutions successives de ces méthodes conduisent à des résultats différents pour 


une même initialisation d'un problème d'optimisation [55]. 


4.4.1. Méthodes déterministes 


Ce sont des méthodes qui s’appuient, pour trouver l’optimum, sur une direction de 
recherche qui peut être fournie par les dérivées de la fonction à optimiser. Elles nécessitent 
donc la connaissance du vecteur gradient. Cette particularité constitue un inconvénient 
majeur dans le cas d’une fonction objectif (multimodal) possédant plusieurs optimums ; 
elles peuvent en effet, converger vers un optimum local. 

La complexité de ces méthodes croit exponentiellement avec la taille de problèmes 
à traiter ; elle devient très importante face à des problèmes à plusieurs variables, fonctions 
objectif et critères. 

Ces méthodes présentent ces inconvénients: 
1. Elles sont souvent inadaptées aux problèmes de grande dimension ; 
2. Elles peuvent être facilement piégées par un optimum local : risque de convergence 
prématurée vers un optimum local. 
3. Elles ne sont pas guidées par un processus intelligent qui conduit la recherche vers 
une bonne solution. 
Il existe de nombreuses méthodes exactes y compris l’algorithme du simplexe, 
méthode du gradient, programmation dynamique, la méthode multistart, la technique de 
programmation quadratique et séquentielle, méthode de Newton ou encore la méthode du 


point intérieur. 


4.4.2. Méthodes stochastiques 


Ce sont des méthodes itératives, directe qui nécessitent n1 point de départ, ni la 
connaissance de gradient de la fonction objectif, pour atteindre la solution optimale. 
Contrairement aux méthodes déterministes, les méthodes stochastiques, possédant un 
comportement aléatoire approximatif, fournissent une solution de qualité satisfaisante en 


un temps de calcul en général raisonnable, sans garantir l’optimalité de la solution obtenue. 
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Elles explorent l’espace des solutions par essais successifs en recherchant les directions les 
plus favorables. Cependant elles demandent un nombre important d’évaluations de la 
fonction objectif en comparaison avec les méthodes déterministes exploitant la dérivée de 
la fonction objectif [28]. 

Parmi les différentes méthodes stochastiques, nous allons uniquement nous 
intéresser aux méta-heuristiques. L’arrivée de cette nouvelle classe de méthodes marque 
une grande révolution dans le domaine de l’optimisation. En fait, elles sont applicables sur 
une grande variété de problèmes d’optimisation de différentes complexités ; à savoir les 
problèmes combinatoires. En outre, Ces méthodes tirent leur intérêt de leur capacité de 
fournir des solutions de très bonne qualité on évitant la convergence prématurée vers un 
optimum local. 

Le mot « heuristique » vient du grec « heurein » « découvrir » et qualifie tout ce qui 
sert à la découverte, à l’invention et à la recherche. Les méta heuristiques progressent de 
façon itérative, où elles manipulent une ou plusieurs solutions à la recherche de l’optimum. 
Les itérations successives doivent permettre de passer d’une solution de mauvaise qualité à 
la solution optimale, en alternant des phases d’intensification, de diversification et 
d’apprentissage. L’état de départ est souvent choisi aléatoirement, ensuite l’algorithme se 
déroule et s’arrête après avoir atteint un critère d’arrêt, généralement en l’atteinte du temps 
d’exécution imparti ou en une précision demandée [50, 59]. 

La majorité des méta-heuristiques sont inspirées des systèmes naturels qui relèvent 
de la physique (recuit simulé qui est inspiré d’un processus métallurgique), de la biologie 
et des principes de l’évolution Darwinienne (les algorithmes évolutionnaires et les 
algorithmes génétiques), de la mémoire des êtres humains (algorithme la recherche tabou), 
ou encore de l’éthologie ou l’intelligence d’essaim (les algorithmes de colonies de fourmis 
ou l’optimisation par essaim particulaire) [52]. La figure 4.3 illustre les différents types 


des algorithmes méta-heuristiques. 
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Fig.4.3. Différents types des algorithmes méta-heuristiques. 


Légende : 

Evolutionary- algorithms: GA: Genetic algorithm; EP: Evolutionary Programming; 
Physics- based algorithms: GSA: Gravitational Search Algorithm 

Swarm-based algorithms: ABC: Artificial Bee Colony; PSO: Particle Swarm Optimizer 


Human- based algorithms: 


4.5. Méthodes d’optimisation proposées. 


4.5.1. Optimisation par essaim de particules (PSO). 


L’optimisation par essaims de particules, ou Particle Swarm Optimazation (PSO) 
en anglais est une technique d’optimisation métaheuristique qui utilise une population de 
solutions candidates pour développer une solution optimale au problème traité. Elle a été 
proposée pour la première fois par Kennedy et Eberhart en 1995 [37]. Cette méthode est 
inspirée du comportement sociologique des animaux qui évoluant en essaims, tels que les 
poissons qui se déplacent en bancs ou les oiseaux migrateurs. En effet, on peut observer 
chez ces animaux des dynamiques de déplacement relativement complexes, alors 
qu'individuellement chaque individu a une intelligence limitée et une connaissance 
seulement locale de sa situation dans l'essaim. Un individu de l'essaim n'a pour 


connaissance que la position et la vitesse de ses plus proches voisins. Chaque individu 
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utilise donc, non seulement, sa propre mémoire, mais aussi l'information locale sur ses plus 
proches voisins pour décider de son propre déplacement. Des règles simples, telles que 
« aller à la même vitesse que les autres », « se déplacer dans la même direction » ou encore 
« rester proche de ses voisins » suffisent à maintenir la cohésion de l'essaim et permettent 
la mise en œuvre de comportements collectifs complexes et adaptatifs. L'intelligence 
globale de l'essaim est donc la conséquence directe des interactions locales entre les 
différentes particules de l'essaim. La performance du système entier est supérieure à la 
somme des performances de ses parties. Kennedy et Eberhart se sont inspirés de ces 
comportements sociaux pour créer l’algorithme PSO [38]. 

La méthode d’optimisation par essaim de particulaire met en jeu une population d’agents 
pour la résolution d’un problème donné. Cette population est appelé essaim. L’essaim est 
composé d’un ensemble d’individus, ces derniers sont appelés particules. Les particules de 
l’essaim représentent des solutions potentielles au problème traité. L’essaim de particules 
survole l'espace de recherche, en quête de l'optimum global. Le déplacement de chaque 


particule est influencé par les trois composantes suivantes. 


1. Une composante physique: la particule tend à suivre sa direction de déplacement 
courante; 

2. Une composante cognitive: la particule tend à se diriger vers le meilleur site par 
lequel elle est déjà passée: 

3. Une composante sociale: la particule tend à se diriger vers le meilleur site déjà 


atteint par ses voisines. 


Dans l’algorithme de PSO, chaque particule est caractérisée par une vitesse 
(velocity) de déplacement et une position (le vecteur solution). Au cours du processus de la 
recherche, chaque particule change sa position dans l’espace de recherche en fonction de sa 
vitesse actuelle, sa position actuelle, sa meilleure position trouvée au cours des itérations 
passées pbest et la meilleure position trouvée par l’essaim gbest. Son déplacement lui 
permet de mettre à jour sa position et sa vitesse de déplacement à chaque itération 
(Fig.4.4). 

L’algorithme PSO comprend plusieurs paramètres de réglage qui permettent d’agir 
sur le compromis exploration — exploitation. L’exploration est la capacité de tester 


différentes régions de l’espace à la recherche de bonnes solutions candidates. 
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L’exploitation est la capacité de concentrer la recherche autour des solutions prometteuses 
afin de s’approcher le plus possible de l’optimum. 


———————ndividual best 


Pbest 


RUOOÏBN _— 


© 
O 4 sno121d 


Previous Velocity 
Fig.4.4. Représentation schématique du déplacement d’une particule. 


Dans un espace de recherche de dimension n, chaque particule i de l’essaim est 
modélisée par son vecteur de position x; = [x:, x, 


déplacement v; = 


.…., Xin] 
= [vis Via … 


et son vecteur vitesse de 
. Vin]. La qualité de sa position est déterminée par la valeur de 


la fonction objective en ce point. Chaque particule garde en mémoire la meilleure position 
par laquelle elle est déjà passée, que l'on note Ppes 


[Pbest1, Dbest2; +, Dies: ainsi que la 
meilleure position atteinte par toutes les particules de l'essaim, notées gpes = [gbesti, best? 
.. Lbestn]- 


En effet, à l’itération t+1, le vecteur vitesse et le vecteur position sont calculés 
respectivement à partir des deux équations suivantes, 


V; (e + 1) =W.V; ()+ C, ®rand, ()® ( pbest;(t)— x;(#) ) (4.9) 
+C; ® rand, (r)@( gbest,(t)— x;(r) ) | 


x;(t+1)= x;(t)+ v, (+1) 


(4.10) 


Le symbole @ signifie ici la multiplication des vecteurs élément par élément, 


LAS L 


pbest, est la meilleure position individuelle de la particule "à 


, gbest; est la meilleure 
position dans le voisinage, V(r) est la vitesse de la particule à à l’itération #, X,(r) est la 


osition de la particule " i "a l’itération f, rand, et rand, sont des nombres aléatoires tirés 
1 > 


79 


uniformément de l’intervalle [0,1] pour chaque itération, C,et C,sont les coefficients 


d’attraction dont les valeurs sont compromises entre 0 et 2, w est le facteur d’inertie. 


e _w.V* correspond à la composante physique du déplacement. Le paramètre w 


contrôle l'influence de la direction de déplacement sur le déplacement futur. 


e CC @rand PE pbest — XF ] correspond à la composante cognitive du 


L 


déplacement. C, contrôle du comportement cognitif de particule. 


e C; @rand, PE gbest* — X} ] correspond à la composante sociale du 


déplacement. C, Contrôle d’aptitude sociale de particule. 


1" 


Le coefficient d’inertie " w "” joue un rôle important dans la convergence de 
l’algorithme. Il est utilisé pour contrôler l’impact que peut avoir la vitesse d’une particule 
sur le comportement de celle-ci, régulant ainsi le compromis entre exploration et 
exploitation. Une valeur élevée de " w " permet une exploration globale de l’espace de 
recherche alors qu’une faible valeur favorise une exploitation locale. De ce fait, il peut être 
envisagé de commencer le processus avec une valeur élevée de w, puis de diminuer "w " 
au fil des itérations, en vue d’affiner la recherche au fur et à mesure de la convergence. 
Cette stratégie a beaucoup amélioré les performances de PSO pour plusieurs problèmes 
d’optimisation. En effet, Le coefficient d’inertie " w " varie linéairement avec le temps 
selon la formule suivante, 


iter 


(4.11) 


max (Wimax Win )- 
LEP ax 


où iter est l’itération courante et iter..,, est le nombre maximal d’itérations. w,,x, et Win 
désignent respectivement les valeurs maximum et minimum du coefficient ” w ". 

L’algorithme d’optimisation par essaim de particules, PSO algorithme, est résumé 
dans les points suivants [40]: 

1. Génération de la population initiale : La population initiale est générée d’une façon 
aléatoire dans l’espace de recherche (dans les limites permises). Chaque solution est 
considérée comme pbest. La meilleure valeur évaluée de pbesrest placée dans 

gbest. 


2. Recherche d’une nouvelle position pour chaque agent : La valeur de la fonction 


objective est calculée pour chaque solution (particule). Si la valeur d’un agent est 
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meilleure que son pbest actuel ; pbest prend cette nouvelle valeur. Si la 
meilleure valeur de pbest est meilleure que gbest actuel ; gbest est remplacé par 
celle-ci et la particule qui correspond à cette valeur est ainsi stockée. 

3. Modification de chaque point de recherche : La solution actuelle est mise à jour en 
utilisant les équations (4.1), (4.2) et (4.3). 

4. Vérification de l’état de sortie: Si le nombre d’itération atteint le nombre 


maximum d’itération iter,,, (ou d’autre critère d’arrêt est satisfait) alors fin du 


programme, autrement, retourner à l’étape 2. 


4.5.2. Méthode de recherche par motifs : Pattern search (PS). 


La méthode de recherche par motifs ou Pattern Search (PS) en anglais est une 
technique évolutionniste à recherche directe (Ne nécessite pas le gradient du problème à 
optimiser) appliquée avec succès pour résoudre une variété de problèmes d'optimisation 
liés à de nombreux domaines de recherche complexes. Cette méthode, introduite par 
Hooke et Jeeves [41], est appartient à la famille des méthodes d’explorations locales c'est- 
à-dire qu’elle consiste, à partir d’une solution initiale " XO ", et par transformations 
successives, à construire une suite de solutions améliorant la fonction objectif. 

Le Pattern Search est caractérisé par une série de déplacements exploratoires autour 
du point (solution) courant formant un maillage (mech). À chaque itération la fonction 
objectif est évaluée sur les points du maillage. Si une amélioration est trouvée, le point 
associé est accepté comme nouveau point courant, et la taille du prochain maillage est 
conservée ou augmentée (Fig.4.5). Sinon, la taille du nouveau maillage, généré autour de 
l’ancien point courant, est réduite. 

Le mécanisme de recherche de l'algorithme PS est représenté par l’organigramme 


présenté sur la figure 4.6. 
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Fig.4.5. Dix premières itérations de l'algorithme Pattern Search. 


Mécanisme de recherche de PS. Les étapes suivantes résument le mécanisme de 


recherche de l'algorithme de Pattern Search [42]: 


1. 


L'algorithme commence par établir un ensemble de points appelés maillage (mech), 
autour des points initiaux " X0 ", ces points initiaux ou points de départ sont soit 
donner par l'utilisateur, soit calculés à partir de l'étape précédente de l'algorithme, 
dans notre cas ces points représentes les résultats obtenue par l’algorithme de PSO. 
A la première itération, avec un scalaire = 1 appelé taille de maillage, les vecteurs 
de motif (pattern) ou de direction sont formés comme suit [0 1], [1 0], [-1 O0] et [0 - 
1]. 

Les points de maillage sont alors construits par l'addition de vecteurs de direction 
aux points initiaux " X0" comme: 

X0 +[0 1], XO + [1 0], XO + [-1 0] et XO + [0 -1] (Fig.4.5). 

L'algorithme calcule la fonction « objectif » à tous les points du maillage jusqu'à ce 
qu'il trouve celui dont la valeur est inférieure à la valeur de la fonction objectif de 


X0. 


4. S'il y a ce point, Alors le sondage (poll) est réussi et le suivant se produit: 


e Le point minimum trouvé devient le point courant X1. 
e La méthode multiplie le maillage actuel par le double. 


e _L'algorithme passe à la prochaine itération. 


82 


5. La valeur X2 sera définie si une valeur de la fonction « objectif » inférieure à celle 
de X1 est trouvée et les itérations continuent. 

6. Si dans une itération particulière, Aucun des points de maillage n'a une valeur de 
fonction objectif inférieure à celle du point initial / courant, dans ce cas le sondage 
est n'est pas réussi et le suivant se produit: 

e  L'algorithme ne change pas le point courant à la prochaine itération. 
e La méthode multiplie le maillage actuel par 1/2. 
e _L'algorithme passe à la prochaine itération. 
7. Les étapes sont répétées jusqu'à ce que l'un des critères d'arrêt suivants soit 
respecté: 
e La taille de maillage est inférieure à la tolérance du maillage. 
e _L'algorithme atteint la valeur maximale des itérations. 
e Le nombre total des fonctions objectifs évaluées par l'algorithme atteint sa 
valeur maximale. 
e La variation de la valeur de la fonction «objectif » est inférieure à la 


tolérance X. 


XO+O I] 


X0+ [-1 0] X0 X0+ [1 0] 


7 ; : x 
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Fig.4.6. Mécanisme de l’Algorithme PS. 


4.5.3. Hybridation d’algorithme PSO avec le PS. 

L’optimisation par essaim de particules (PSO) est une méta-heuristique qui permet 
de mieux appréhender les grands espaces de recherche mais elle n’a pas un fort pouvoir 
d’intensification. Inversement, la méthode de PS a une très grande capacité 
d’intensification mais cela peut entraîner le blocage de la méthode dans des optima locaux. 
Par conséquent la coopération de ces deux méthodes permet de combiner les avantages des 
deux algorithmes et de maintenir un bon équilibre entre l’exploitation des solutions et 
l’exploration de l’espace de recherche. Donc L'objectif principal de la technique proposée 
est d'améliorer la propriété de recherche globale de PSO par une recherche locale de 
robustesse de PS. L'algorithme hybride s'exécute en deux phases. Dans la première phase, 
l’algorithme PSO commence par l’initialisation de la taille de l’essaim et des positions (x; 
= Xi], Xp, …, Xn) aléatoires des particules de l’essaim. Ensuite, il initialise la meilleure 
position rencontrée par chaque particule (pes = Pbesti> Pbest> +, Pbesm) Par Sa position 


initiale x;. Puis, en se basant sur la fonction objectif, l’algorithme calcule la qualité (fitness) 
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de chaque particule pour pouvoir tirer la meilleure position (solution) rencontrée par 
l’essaim grey. Après chaque génération, l’algorithme PSO envoie chaque solution au 
Pattern Search (deuxième phase) afin d’être améliorée. Les meilleures solutions obtenues à 
partir de PS seront les solutions de PSO pour la prochaine génération. L'échange entre PSO 
et PS se poursuit jusqu'à ce que le critère d’arrêt soit satisfait. La solution finale de la 
première méthode devient la solution initiale de la méthode suivante. Cette enchaînement 
de méthodes peut notamment permettre de procéder à une phase de diversification pour 
couvrir une grande partie de l’espace de recherche (PSO) avant d’arriver à une phase 
d’intensification pour améliorer les solutions obtenues (PS). L'organigramme de l’hybride 


PSO-PS Algorithme est présenté dans la figure 4.7. 


Read Power System Data 
Implémentation of FACTS Devices 
Run Power Flow 


Convergence 


Yes 
Set PSO Result as Starting 
Point XO0 for PS 


Convergence 


Yes 
Save the Final Result 


Fig.4.7. Algorithme hybride PSO- PS. 
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4.5.4. Algorithme Jaya. 


Jaya est un nouvel algorithme d'optimisation stochastique à base de population de 
solutions, développé par Rao en 2016 [43], pour résoudre des problèmes d'optimisation 
avec et sans contraintes. Cet algorithme est inspiré par le fait que les candidats de la 
population devraient se déplacer vers la meilleure solution de la population et éviter la 
mauvaise (le pire). 

L'algorithme tente toujours de se rapprocher du succès (c.-à-d. d'atteindre la 
meilleure solution) et tente d'éviter l'échec (c.-à-d. s'éloigner de la mauvaise solution). 
L'algorithme s'efforce de devenir victorieux en atteignant la meilleure solution et, par 
conséquent, il est nommé Jaya (un mot sanskrit signifiant victoire) [44]. 

Soit " f(x) " la fonction objective à minimiser. À chaque itération " i ", supposons 


U " 


qu'il existe " m nombre de variables (j = 1,2, …, m) et "n 


" 


nombre de solutions 
candidates (taille de la population, k = 1,2, …, n). Soit le meilleur candidat " best" obtient 


la meilleure valeur de F (x) (c.-à-d. f(x), ) dans toutes les solutions candidates " n "et le 
mauvaise candidat " worst " obtient la mauvaise valeur de f (x) (c.-à-d. f(x),,) dans 
l'ensemble des solutions candidates. 

Si X ;,, est la valeur de la j-ième variable pour le k-ème candidat au cours de la i- 
ème itération, alors cette valeur est modifiée selon l'équation suivante [44, 45]: 
bé —|X ju ] (4.12) 


! — 
Lei =X ji + A x 


j,best,i jk,i 


Où x . est la valeur de la variable " j " pour le meilleur candidat et X 


j,best,i j,worst i 


est la valeur de la variable ” j " pour le mauvaise candidat. X°,; est la valeur mise à jour 


de X 4: ñ,;; et ,,,s0nt les deux nombres aléatoires pour la j-ème variable pendant la i- 


| 


indique la tendance de la solution à se rapprocher de la meilleure solution et le terme 


ie 


ème itération tirés uniformément de l’intervalle [0, 1]. Le terme vx res be 


-|x 


) indique la tendance de la solution à éviter la mauvaise solution. 


L me es — IX ei 
Les nombres aléatoires ñ;; et r,;; assurent une bonne exploration de l'espace de 


recherche. La valeur absolue de la solution candidate IX améliore encore la capacité 


j ki 
d'exploration de l'algorithme. 

L'objectif principal de l'algorithme de Jaya est d'améliorer la valeur de ” f{x) " de 
chaque solution candidate dans la population. Par conséquent, l'algorithme tente de 


déplacer la valeur de la fonction objectif de chaque solution vers la meilleure solution en 
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mettant à jour les valeurs des variables. Une fois les valeurs des variables mises à jour, les 
solutions mises à jour sont comparées aux anciennes solutions correspondantes et seules 
les bonnes solutions (c.-à-d. les solutions ayant une meilleure valeur de fonction objective) 
sont considérées pour la prochaine génération. 

Les étapes principales de processus de recherche de l’algorithme Jaya sont 


résumées dans l’organigramme présenté dans la figure 4.7 et expliquées dans ce qui suit. 


Étape 1: Initialisation: Dans cette étape, une population de " n "candidats ayant une taille 
prédéfinie est générée au hasard dans l'espace de recherche. 

Étape 2: Identification des meilleures et des mauvaises solutions: Après l'initialisation de 
la population, les meilleurs et les mauvais candidats parmi la population en 
termes de fonctions objectives sont identifiés. 

Etape 3: Dans cette étape, les candidats de la solution sont déplacés selon la formule 4. 
Pour un candidat donné, si la nouvelle solution (après le déplacement) est 
meilleure que l'ancienne (avant le déplacement), la solution est acceptée et le 
candidat correspondant est mis à jour. Sinon, la nouvelle solution est rejetée et 


l'ancienne est conservée. 
Le processus de l'étape 2 et de l'étape 3 est répété jusqu'à ce qu'un critère d’arrêt 


soit satisfait. Le critère d’arrêt choisi dans ce travail est le nombre maximal d'itérations et il 


est clair que ce critère peut être changé et n'importe quel autre critère peut être implémenté. 
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Initialization 


Identify the best and the worst solutions in 
the population 


Update the solutions of candidates based 
on the best and The worst solutions 
according to equation (20) 


No Yes 


If the new solution for a condidate 
is better than the old or 


Keep the old 
solution 


Accepte the new solution and replace 
the old one in the population 


No Yes 


If the termination criterions 1s 
satisfied 


Report the optimal solution found 


Fig.4.8. Diagramme de l'algorithme Jaya. 
Traduction: 
- Initialisation 
- Identifier les meilleures et les mauvaises solutions dans la population 
- Mettre à jour les solutions des candidats basés sur le meilleur et la mauvaise solution 


selon l'équation (4.12) 


! = 
7m ui La j,ki 


x ) (4.13) 


le M, ji bone _ IX ei 
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- Si la nouvelle solution pour un candidat est meilleure que l'ancienne 
- oui : Accepter la nouvelle solution et remplacer l'ancienne dans la population 
- non : Garder l'ancienne solution 

- Si les critères d’arrêt sont satisfaits 


- oui : exposé la solution optimale trouvée 


L’algorithme développé : Multi-objectifs Jaya (MOJaya). Les principales étapes 
de l'algorithme MOJaya sont: 


Etape 1: Initialisation: dans cette étape, une population initiale notée Poy(0) est générée 
aléatoirement dans l'espace de recherche. De plus, une archive externe vide Popyrchive(O) 
est créé. 
Etape 2: Affectation de la fitness: les valeurs de la fitness pour tous les candidats d'une 
population globale (total) Pop= PopU PopyrchiveCOmposé de Popüir) et PoParchive@t) Sont 
calculées en utilisant la procédure suivante: 

Pour chaque candidat «1», dans la population Popif et dans les archives 
PoOPyrchiveGD une force S(i) est assignée, S(i) représente le nombre de candidats qu'elle 
domine. Sur la base de Si), la fitness brute (the raw fitness) R() d'un candidat «1 » est 


calculée comme suit: 


R(i) = DSC (4.14) 
jePop 


j dominei 
Il convient de mentionner que, puisque le problème de l’écoulement de puissance 
optimal multi-objectifs est formulé comme un problème de minimisation, R(i)=0 
correspond à un candidat non-dominé (une solution du front de Pareto) tandis qu'une 
valeur élevée de R(i) est synonyme que ce candidat est dominé par de nombreux candidats. 


Après cela, une densité D(i) est calculée en utilisant l'expression suivante: 


D(i) = 


(4.15) 


GO; +2 
Où of est le k°7° plus proche voisin (basé sur la distance dans l’espace des objectifs), 


k= nn .Le 2 est ajouté dans l'expression de D(i) pour s'assurer que celui-ci est compris 


entre 0 et 1. 
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Finalement, la fitness F(i) est calculé en utilisant l'expression suivante: 
F(i) = R()+ D() (4.16) 


Étape 3: Sélection environnementale : Dans cette étape, les candidats non-dominés dans 


Pop) et dans l'archive PopyrchiveGf) Sont copiés dans Popyrehive@t+1). S1 la taille de 
POParchiveGt+1) est supérieure à « na », POPurchive@t +1) est réduit en utilisant l'opérateur 
de troncature (the truncation operator), sinon, si la taille de Popyrchive@t+1) est inférieure à 
«na», alors PoPyrchive(@it+1) est complété par les candidats moins dominés dans Pop{ir) 
et dans PoPyrchive@n) . Cela permettra à la taille de Popyrchive d'être constante sur toutes les 


itérations quel que soit le nombre de candidats non-dominé trouvés. 


Étape 4: Identification des meilleures et des mauvaises solutions dans Pop basé sur la 

fitness F. 

Étape 5: Déplacer les candidats de Pop en utilisant l’équation (4.14). Ensuite, pour 

chaque candidat dans Pop, si une nouvelle solution (après déplacement) domine l'ancienne 

(avant déplacement) cette solution est mises à jour. Sinon, la nouvelle solution est rejetée. 
Ce processus de l'étape 2 à l'étape 5 est répété jusqu'à ce qu'un critère de terminaison 


soit atteint. 


4.6. Technique de la logique floue pour l'extraction de la solution de compromis. 


Une fois que l'ensemble optimal de Pareto est obtenu, le décideur doit choisir une 
solution de compromis parmi les solutions optimales qui satisfait leur cahier de charge. 
Cette opération reste très complexe. En raison de la nature imprécise du jugement du 
décideur. Pour la prise de décision, c'est-à-dire pour aider le décideur, il est nécessaire 
d'identifier la meilleure solution de compromis parmi l'ensemble de solutions non-dominée 
trouvées. Dans ce travail de thèse, nous proposons d’utilisation de la théorie logique floue 
(Fuzzy logic) pour la sélection de la meilleure solution de compromis déterminée comme 
suit. Tout d'abord, chaque fonction objectif de la solution non-dominée est remplacée par 
une fonction d'appartenance. Les fonctions d'appartenances présentent les buts de chaque 
fonction objectif ; dont elles sont définies par les expériences et la connaissance intuitive 
du décideur. Dans ce cas, nous avons choisi une fonction d'appartenance linéaire simple 


définie comme (Fig.4.9) [36, 40]: 
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Es Fax 


Fig.4.9. Fonction linéaire d’appartenance de la fonction objectif. 


1 F; FN da 
es pin 
L 


max min 
F; — F; 


Uj= SRSFT (4.17) 


0 FAR ep 
avec FMitet FM&sont la valeur minimale et maximale de la i°"* fonction objectif parmi 
toutes les solutions non dominées, respectivement. La fonction d’appartenance « u » est 
variée entre le 0 et le 1, avec H=0 indique l'incompatibilité de la solution avec l'ensemble, 
tandis que H=1 signifie la pleine compatibilité. 
Ensuite, pour chaque solution non-dominée « k », la fonction d'appartenance normalisée 


ul est calculée comme suit: 


Nobj k 
K 2 max Di É 4.18) 
# D Nobj k ( j 
ei Zi Fi 


avec N,5j €t Nas Sont le nombre des fonctions objectives traitées et le nombre des 
solutions non-dominées trouvées, respectivement. Donc, la meilleure solution de 


compromis est celle qui à la grande valeur de u* . 
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CHAPITRE 5 . SIMULATIONS ET RESULTATS. 


5.1. Réseau test IEEE 30-bus. 


Le réseau électrique test IEEE 30-bus est utilisé dans la littérature pour évaluer des 
méthodes ou des modèles développés dans des divers travaux de recherche. Ce réseau, 
choisi pour notre simulation, est un réseau standard qui représente une portion du système 
de puissance électrique Américain. Le schéma unifilaire correspondant est représenté sur la 
figure 5.1. Il est constitué de trente (30) jeux de barres et six (6) générateurs connectés aux 
nœuds 1, 2, 5, 8, et 13. Ces génératrices alimentent vingt (20) charges (PQ bus) à travers 
quarante et une (41) lignes de transport, et quatre (4) transformateurs à prises de réglage, 
reliant, à travers des lignes, les jeux de barres (4-12), (6-9), (6-10) et (28-27). En plus, deux 
batteries de condensateurs fixes sont installées aux jeux de barres 10 et 24. La puissance 
totale demandée, respectivement active et réactive, est de 283.4 MW et 126,2 MVAR. Les 
valeurs de tension et de puissance de base sont U,= 135 kV et S,=100 kVA respectivement. 
Les données détaillées du réseau test IEEE 30-nœuds sont donnés dans les fableaux S5.I et 
5.2 [60]. Les coefficients du coût et les limites inférieures et supérieures des puissances 


actives et réactives des générateurs sont données dans le fableau 5.3, [61]. 


Fig. 5.1. Schéma unifilaire du réseau électrique IEEE 30-bus. 
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Tab.5.1. Données nodales du réseau (30 nœuds). 


Code | Mag. | Deg. a Load No | Code | Mag. | Deg. nine Load 
RE Rat uw La 


© ® I A 1 BR © D = 


= En mi it Hi pt et ed ot et pt et et pt ji 


0,0264 
0,0204 
0,0184 
0,0042 
0,0209 
0,0187 
0,0045 
0,0102 
0,0085 
0,0045 


© 


S.0dOHHOTTSOOOOCC 


[SMW)hr] 
0.00375 
0.01750 
0.06250 
0.008334 
0.02500 
0.02500 


= 

LEA 
4 

Led 


DOSCOCCSe 
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Tab. 5.4. Valeurs limites des variables de contrôle 


Variables 
Tensions des générateurs VG (p.u.) 
Rapports des transformateurs T (p.u.) 


Puissances réactives des compensateurs shunts | Qc (p.u.) 
Réactance du compensateur série Xrcsc (p.u.) 
Susceptance du compensateur shunt Qsvc (MVAr) 


Les valeurs limites des modules de tensions des générateurs (PV-bus), des rapports 
des transformateurs, des puissances réactives des compensateurs shunts, de la réactance du 
compensateur série TCSC et de la susceptance du compensateur shunt SVC sont données 
dans le tableau 5.4. Les limites inférieures et supérieures des modules de la tension de tous 
les jeux de barres de charge (PQ-bus) sont choisies comme  [0.95; 1.1 p.u.], 


respectivement [61]. 


5.2. Optimisation Multi objectif de l’écoulement de puissance par L'algorithme 
hybride PSO-PS avec intégration des dispositifs FACTS. 


Dans cette partie de travail, on considère cinq fonctions « objectifs »; la 
minimisation du coût de combustible, l’amélioration du profil de tension, l’amélioration de 
la stabilité de tension, la minimisation de la déviation de tension et la minimisation des 
pertes actives et réactives de transmission en présence du système FACTS. Parmi les 
dispositifs FACTS, on choisi ceux qui sont plus conformes aux applications spécifiques 
considérés ; tels le maintien de la tension et le contrôle de l’écoulement de puissance. Les 
dispositifs FACTS considérés dans ce travail sont le SVC et le TCSC. Le SVC étant le plus 
utilisé, puisque performant, dans la compensation de l’énergie réactive et le maintien de la 
tension ; le TCSC représente le meilleur choix quant au contrôle approprié et à la bonne 
gestion de l'écoulement de puissance et, par suite, à la réduction des pertes de puissance 
dans le réseau [31]. L'emplacement optimal des SVC et TCSC est pris de la référence [62] 


et est spécifié dans le tableau 5.5. 


Tab.5.5. FACTS devices location. 


Dispositifs FACTS | Emplacement optimal 


TCSC Branche: 28-27 (Line 36) 
SVC Noeud: 22 
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Dans ce qui suit, on traite, particulièrement, de l'impact des dispositifs FACTS sur 
les différents types des fonctions « objectifs » retenues, à cet effet, dans ce travail. 
Les paramètres de l'algorithme PSO-PS pour l'optimisation sont donnés dans le tableau 


5.6. L'optimum global est le meilleur résultat après 10 exécutions du programme. 


Tb. 5.6 Simulation parameters 


parameters | PS parameters 


Population size mesh size 1 
Number of generation Number of generation 150 


Wimax / mesh expansion factor 2 
Win : mesh contraction factor 0.5 
C1, C2 


CAS. 1 : Minimisation du coût de combustible. 


L'objectif dans ce cas est de minimiser le coût de combustible pour tous les 
générateurs. Les variables optimales de contrôle, obtenues à partir de PSO et PSO-PS 
hybride en absence et en présence de dispositifs FACTS sont présentées dans le tableau 
5.7. Ce dernier présente les résultats concernant le coût minimal de production, la déviation 
de la tension, l'indice de stabilité et les pertes de puissance. Ces résultats montrent que le 
coût minimum obtenu est amélioré de 799.9738 $ / h, en utilisant PSO, à 7998723 $ /h, 
en utilisant PSO-PS mais le meilleur coût, celui de 799.5485, est réalisé suite à l'insertion 
des dispositifs FACTS. Ici, 1l est important de noter que les dispositifs FACTS peuvent 
réduire le coût total du combustible dans le système électrique. A ce coût réduit minimal, 
en présence des FACTS, correspond une déviation de tension de 1,4590 p.u., des pertes 
de puissance active et réactive totales ainsi qu’un indice de stabilité de la tension de 8,7618 
MW, - 2,169 MVAR et 0,1214 p.u., respectivement. Les caractéristiques de convergence, 
déduites de la minimisation du coût de combustible, ainsi que les approches proposées, 
sont illustrées en fonction des itérations sur la figure 5.2. En outre, les résultats obtenus 
sont comparés à ceux d'autres méthodes appliquées dans divers travaux spécifiques à la 
problématique posée, (Tb.5.8). On en déduit, clairement, que la méthode PSO-PS proposé 


donne le meilleur minimum global par rapports aux autres méthodes. 
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Tb.5.7. Optimal control variable settings for minimization of fuel cost case 
178.1960 | 177.1504 177.0989 
48.9732 | 48.7208 48.7088 
21.390932 | 21.3146 21.3093 
21.7776 | 21.1636 21.1228 
Peu 10.0000 | 11.9056 11.9211 
Pc13 12.0000 | 12.0000 12.0000 
Va pu.) 1.1000 1.1000 1.1000 
1.0879 1.0879 1.0878 
1.0615 1.0616 1.0616 
1.0696 1.0697 1.0694 
Vou 1.1000 1.1000 1.1000 
Vois 1.1000 1.1000 1.1000 
T1 (6.9) 1.0468 1.0465 1.0668 
T2 (610) 0.9000 0.9000 0.9000 
T5 (4-12) 1.0348 1.0343 1.0024 
T6 28-27) 0.9747 0.9746 0.9674 
QCio (ù 0 
QC 
QCis 
QCr 
QC 
QCu 
QC: 
QC 
QC 
X_TCSC 0.3131 
Q_SVC 18.1168 
Cost ($/hr) (*) 799.5485 
VD (p.u.) 1.4590 
L_max (p.u.) 0.1214 
Ploss MW) 8.7618 
Qloss (MVar) -2.1629 


(*) : Meilleures valeurs de la fonction « objectif » (en gras). 


So © © © © © © © 
S © © © © © © © 
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Tb.5.8. Comparison of results for minimization of fuel cost case 


Cost (Shr) 


PSO-PS with 


FACTS 799.5485 - 


EGA-DQLF 799.56 (Kumari and Maheswarapu [63] 
799.8723 - 
799.9217 (Bouchekara 2014) [36] 
800.078 (Bouchekara 2016) [64] 


800.0831 (Surender Reddy et al. 2014) [65] 
800.3887 (Abaci and Yamacli 2016) [66] 
800.9617 (Mahdad et al. 2010) [67] 
801.5733 (Duman 2016) [68] 
801.75 (Daryani et al. 2016) [69] 
801.81 (Roselyn et al. 2011) [70] 
801.96 (Niknam et al. 2012b) [71] 


807 : : : - 

806 He PSO-PS | 
CELL LE PSO -PS “| 

805 with FACTS 


Cost ($/hr) 


0 100 200 300 400 500 
lterations 


Fig. 5.2. Fuel cost variations with different proposed approaches. 
CAS. 2 : Minimisation du coût de combustible et amélioration du profil de tension. 


La tension du jeu de barres est un indice fondamental déterminant de la qualité et 
de la sécurité dans les systèmes d'énergie électrique. Minimiser seulement le coût de 
combustible peut impliquer des états de profil non admissibles de la tension. Par 
conséquent, une fonction supplémentaire d’amélioration du profil de tension est introduite 
en somme avec la fonction de coût du combustible. L'objectif est de réduire le coût du 
combustible et d'améliorer le profil de tension simultanément en maintenant les valeurs de 
tension du jeu de barres PQ proches de 1,0 par unité. La fonction « objectif » peut être 


exprimée comme [61]: 
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F = Ecost + W. Voltage __deviation (5.1) 


où, w est un facteur de pondération approprié; Fos and Fvyonage deviarion SON SpéCifiées 


précédemment comme suit, 


NG 
RE 1; (5.2) 
i=1 
FVotage _ deviation — > 4 Il (5.3) 
ieNL 


Les résultats obtenus en utilisant PSO-PS avec et sans dispositifs FACTS sont 
saisies dans le tableau 5.9 et comparés à ceux d'autres méthodes ; telles que BHBO et EM, 
dans le tableau 5.10. La courbe de convergence des résultats de déviation de tension en 
fonction des itérations pour différentes approches est illustrée sur la figure 5.3. La 
variation du coût du combustible et de la déviation de tension pour l'algorithme PSO-PS en 
présence de dispositifs FACTS est illustrée sur la figure 5.4. On peut voir d'après le 
tableau 5.9 que la déviation de tension optimale obtenue est améliorée de 0,1043 p.u., en 
utilisant PSO, à 0,0984 p.u., en utilisant PSO-PS et réduite à 0,0955 p.u. après insertion 
des dispositifs FACTS. Il convient, ainsi, de mentionner dans ce cas que les dispositifs 
FACTS sont très efficaces pour la minimisation de la déviation de tension. En outre, la 
déviation de tension obtenue dans ce deuxième cas est fortement réduite par rapport celle 
du cas précédent, (réduite de 1,4590 p.u. à 0,0955 p.u. ; soit de % 93,45). Les amplitudes 
de toutes les tensions nodales, étant en même temps, dans les limites requises de sécurité. 
Le coût obtenu en présence de FACTS pour ce cas est de 804.3092, tandis que les valeurs 
de Liax Ploss €t Qioss Sont de 0.1451 p.u., de 99887 MW et de 7.7490 MVAR, 
respectivement. De plus, la valeur de la puissance réactive optimisée délivrée par le SVC 
est de 10,2256 MVAR, tandis que celle de Xrcsc est de 0,6036 p.u. D'après les résultats 
obtenus  (Tb.5.10), il est clair que l'approche proposée donne de meilleurs résultats 


comparativement à d'autres méthodes. 
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Tb. 5.9. Optimal control variable settings for voltage profile improvement case 


Control Variable PSO-PS | PSO-PS with FACTS 


Pa MW) 184.8192 | 176.5139 176.0196 


Po 50.9818 48.8777 48.8563 
Pcs 22.3303 21.8114 21.6847 


Ps 10.0000 | 22.2434 22.5944 
Peu 13.9845 | 11.8459 12.2361 
Pc 12.0000 |  12.0000 12.0000 
Vai (pu) 1.0414 1.0365 1.0320 
Vo 1.0243 1.0214 1.0153 
1.0153 1.0116 1.0146 
Vos 1.0013 0.9996 1.0001 
Veu 1.0758 1.0789 1.0532 
Vas 0.9865 0.9938 1.0070 
Les 1.1000 1.1000 1.0689 
Das 0.9000 0.9000 0.9108 
Lsdis 0.9422 0.9535 0.9717 
Éécéen 0.9475 0.9654 0.9661 
QCi0 3.6696 0 
QC» 2.4767 0.3305 
QCis 5.0000 4.9998 
QC: 0 0.9336 
QC 5.0000 5.0000 
QCA 5.0000 0.1169 
QC: 4.9991 4.7983 
QC 5.0000 5.0000 
QC 2.2510 3.0904 
X_TCSC 2 0.6036 
Q_SVC - - 10.2256 
Cost ($/hr) (*) | 805.7064 | 803.8194 804.3092 
VD(pu) (* 0.1043 0.0984 0.0955 
L_max (p.u.) 0.1380 0.1369 0.1451 
Ploss (MW) 10.7132 |  9.8899 9.9887 
Qloss (MVar) 11.2352 |  8.7207 7.1490 


(*) : Meilleures valeurs de la fonction « objectif » (en gras) 


Tb. 5.10. Comparisons of results for voltage profile improvement case 


Voltage Deviation (p.u.) 


PSO-PS with FACTS - 


PSO-PS - 
BHBO (Bouchekara 2014) [36] 
EM (Bouchekara 2016) [64] 


99 


© 
[ep] 
L 


PSO 
a  — PSO -PS - 
CS PSO - PS 
_. with FACTS 
5 0.4 
C 
2 
0 
(o) 
TD 
(o] 
g 02 
o) 
> 
0.1 
0 r r L L Æ 
(0) 100 200 300 400 500 
lterations 
Fig. 5.3. Voltage deviation minimization with 
different proposed approaches. 
880 É c c c 0.8 
ELLE CELL Cost 
+ Voltage deviations || — 
= 
a 
[®] 
ue C 
Eu © 
& E 
” ë 
8 d 
O (o) 
D 
œ 
+ ER IEEE TR > 
_ ne 0 nn don om Qu me mm men mn 1 
CRETE] ass DLL OPPOSER EDP ELEL LESC EL CELL IS CECELELEES LEPELEER 
800 7 - - | #0 
0 100 200 300 400 500 


lterations 


Fig. 5.4. Variations of fuel cost and voltage deviations for case 2 


CAS. 3 : Minimisation du coût de combustible et amélioration de la stabilité 


de tension. 


Dans ce cas, l'objectif est de minimiser le coût du carburant et d'améliorer la 
stabilité de la tension simultanément. La fonction objective choisie peut être formulée, 


alors, par l'équation suivante: 
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F=F cost + w.F Voltage _ stability _enhancemern (54) 
avec, 
NG 
Ecost De (5.5) 
i=1 
F Voltage _ stability __enhancement — Lnax (5.6) 


où, w est un facteur de pondération approprié. 


Les valeurs optimales des variables de contrôle et des fonctions « objectifs » 
relatives à l'amélioration de la stabilité de tension sont données dans le Tableau 5.11. Ici, 
il est important de confirmer que toutes les contraintes de sécurité telles que les tensions 
aux nœuds de charge et les puissances de lignes sont satisfaites et sont dans les limites 
admissibles. L’évaluation comparative des résultats, concernant la valeur de l'indice de 
stabilité, obtenus par les approches proposées au cours du processus de recherche est 
illustrée sur la figure 5.5. Une comparaison de ces résultats avec différentes méthodes est 
présentée dans le tableau 5.12. La courbe de convergence du coût de combustible et de 
l'indice de stabilité pour l'algorithme PSO-PS en présence de dispositifs FACTS est 
présentée sur la figure 5.6. Les résultats du tableau 5.11 montrent que le coût total du 
combustible et l'indice de stabilité Lmax déterminés par l'algorithme PSO-PS en présence 
de FACTS sont respectivement de 799,9284 $ / h et de 0,0758 p.u. ; ce qui montre, donc, 
une réduction de 37,56% du L,., et une augmentation d'environ 0,047 % du coût total 
du combustible par rapport au cas 1 (coût: 799,5485 $ / h, Lmax: 0,1214 p.u.). En 
comparant les résultats obtenus, respectivement, en absence et en présence de dispositifs 
FACTS, on voit que l'introduction de dispositifs FACTS peut considérablement améliorer 
l'indice de stabilité (33.97%), (de 0.1148 p.u à 0.0758 p.u.). De plus, ces résultats sont 
significativement diminués par rapport au cas 1 (0,1214 p.u.) et au cas 2 (0,1451 p.u.). De 
même, la puissance réactive optimisée délivrée par le SVC, pour ce cas, est de 22.4861 
MVAR, tandis que les pertes de puissance active et réactive sont 8.7801 MW et -3.9749 
MVar, respectivement. Le tableau 5.12 montre que le résultat obtenu est meilleur par 


rapport à ceux autres méthodes. 
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Tb. 5.11. Optimal control variable settings for Voltage stability enhancement case. 


12.0000 12.0000 12.0000 


1.0834 1.1000 
1.0888 


| 10046 | 1.0300 1.0725 


1.1000 1.1000 


1.1000 1.1000 
| 09000 | 0.9000 1.0676 


5.0000 


5.0000 
5.0000 


5.0000 


Q_SVC 


= 


(#) : Meilleures valeurs de la fonction « objectif » (en gras). 


102 


Tb. 5.12. Comparison of results for voltage stability enhancement case 


0.10402 (Kumari and Maheswarapu 2010) [63] 


0.1052 | (Surender Reddy et al. 2014) [65] 
PSO-PS DST 
0.1156 | (Bouchekara 2016) [64] 
BHBO 0.11671 | (Bouchekara 2014) [36] 
0.1244 | (Abaci and Yamacli 2016) [66] 
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Fig.5.5. Stability index minimization with 
different proposed approaches. 
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Fig. 5.6. Variations of fuel cost and stability index for Case 3. 
CAS. 4: Minimisation des pertes actives. 


Les paramètres optimaux des variables de contrôle et les valeurs des fonctions 


objectives obtenus par l’algorithme PSO, PSO-PS en l'absence et en présence de dispositifs 
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FACTS sont montrés dans le fableau 5.13; sa caractéristique de convergence est donnée 
sur la figure 5.7. Il est important de confirmer, ici également, que toutes les variables de 
contrôle et d'état sont dans leurs limites requises. D'après le fableau 5.13, la valeur 
optimale des pertes actives obtenue par l'algorithme PSO-PS en présence du dispositif 
FACTS est de 2.7914 MW, laquelle est réduite d'environ 2,097% par rapport aux 2,88512 
MW obtenus par PSO-PS (sans FACTS) et réduite d'environ 12,08% par rapport aux 
3,1752 MW obtenus par PSO; de plus cette valeur est meilleure par rapport à celle 
obtenue par différentes autres méthodes (Tb.5.14). Dans ce cas, la puissance réactive 
optimisée délivrée par le SVC est de 5,99698 MVAR, tandis que Xrcsc; Lmax €t Qioss 
sont de -0,2487, de 0,0827 p.u. et de -24.7175 MVAR, respectivement. 


Tb. 5.13. Optimal control variable settings for real power losses case 


Control Variable PSO-PS | PSO-PS with FACTS 


51.5757 | 51.2498 51.1920 
80.0000 | 80.0000 80.0000 
50.0000 | 50.0000 50.0000 
35.0000 | 35.0000 35.0000 
Peu 30.0000 | 30.0000 30.0000 
Pc 40.0000 | 40.0000 40.0000 
Va (pu) 1.0331 1.1000 1.1000 
1.0281 1.0976 1.0976 
1.0067 1.0799 1.0799 
1.0129 1.0869 1.0868 
Vu 1.1000 1.1000 1.1000 
Vous 1.1000 1.1000 1.1000 
T1 (6.9) 0.9000 1.0550 1.0690 
T12 (6-10) 0.9000 0.9000 0.9000 
T5 412) 0.9331 0.9842 0.9841 
T6 28-27) 0.9136 0.9727 0.9824 
QC: 0 5.0000 5.0000 
QC2 5.0000 5.0000 5.0000 
QCi15 5.0000 5.0000 5.0000 
QC: 5.0000 5.0000 5.0000 
QC 0 5.0000 4.1859 
QC 5.0000 5.0000 5.0000 
QC: 5.0000 3.5972 2.7385 
QC 5.0000 5.0000 5.0000 
QC 2.8784 2.5213 2.4170 
X_TCSC - -0.2487 
Q_SVC - 5.9698 
VD (p.u.) 2.0503 2.0825 
L_max (p.u.) 0.1150 0.0827 
Ploss MW) (*) 2.8512 2.7914 
Qloss (MVar) -24.3696 -24.7175 
(*) : Meilleures valeurs de la fonction « objectif » (en gras). 
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Tb.5.14. Comparison of results for real power loss minimization 


SOS 2 
6.5 t t r = 
PSO 
———— PSO - PS 
5.5 CÉLLLLLELEN . .. ee 


PLosses (MW) 


r Fa E DT 
0 100 200 300 400 500 
Iterations 


Fig.5.7. Real power losses minimization with 
different proposed approaches. 


CAS. 5: Minimisation des pertes réactives. 


Les valeurs optimales des variables de contrôle et des fonctions objectives sont 
donnés dans le fableau 5.15 On en déduit, de ce dernier, que la valeur la plus optimale des 
pertes réactives, obtenues, alors, par l'algorithme PSO-PS en présence du dispositif 
FACTS est de -26.0946 MVAR ; laquelle valeur est réduite d'environ 3,41 % par rapport 
à celle (-25,2040 MV AR) obtenue par PSO-PS sans FACTS et d'environ 7,45% par rapport 
à celle (-24,1499 MV AR) obtenue par PSO et est inférieure aux valeurs obtenues dans tous 
les cas précédents. Ces résultats sont meilleurs en comparaison avec les résultats obtenus 
avec les méthodes EM et BHBO, comme le montre le tableau 5.16. De plus, la puissance 
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réactive optimisée délivrée par le SVC dans ce cas est de 9,8566 MVAR, tandis que Lx et 
Pioss Sont de 0,0770 p.u. et 2.8590 MW, respectivement. La différence d’effet sur les pertes 
réactives, impliquée par l’utilisation des PSO, PSO-PS avec et sans FACTS, est donnée sur 


la figure 5.8. 


Tb.5.15. Optimal control variable settings for reactive power losses case 


Control Variable PSO-PS | PSO-PS with FACTS 


51.4581 | 51.3087 51.2577 
80.0000 | 80.0000 80.0000 
50.0000 | 50.0000 50.0000 
35.0000 | 35.0000 35.0000 
Pau 30.0000 | 30.0000 30.0000 
Peu 40.0000 | 40.0000 40.0000 
Va (pu) 1.1000 1.1000 1.1000 
Va 1.1000 1.1000 1.1000 
Ves 1.0920 1.0920 1.0921 
Ves 1.1000 1.1000 1.1000 
Vau 1.1000 1.1000 1.1000 
Vas 1.1000 1.1000 1.1000 
T1 (6-9) 0.9650 1.0016 1.0004 
T2 (610) 1.1000 0.9660 1.0034 
T5 4.12) 0.9757 0.9948 0.9941 
T6 28-27 0.9796 0.9861 0.9983 
QCi10 0 5.0000 5.0000 
QC 5.0000 5.0000 5.0000 
QCi5 0 5.0000 5.0000 
QC:r7 5.0000 5.0000 5.0000 
QC 5.0000 5.0000 5.0000 
QC: 0 5.0000 5.0000 
QC: 5.0000 5.0000 3.9636 
QC 5.0000 5.0000 5.0000 
QC» 4.0341 3.2714 3.4325 
X_TCSC - - 0 
Q_SVC - - 9.8566 
VD (p.u.) 1.8736 2.1318 2.1679 
L_max (p.u.) 0.117 0.1142 0.0770 
Ploss (MW) 3.0571 2.9102 2.8590 
Qloss (MVar) (*) | -24.1499 | -25.2040 -26.0946 
(*) : Meilleures valeurs de la fonction « objectif » (en gras). 
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Fig.5.8. Reactive power losses minimization 
with different proposed approaches. 


Tb.5.16. Comparison of results for reactive power losses minimization 


Method Reactive Power Loss MVAR) 


PSO-PS with FACTS -26.0946 - 


-25.204 - 
—22.0196 (Bouchekara 2016) 
—20.1522 (Bouchekara 2014) 


Commentaire. 


Dans cette première partie, l'algorithme d'optimisation hybride PSO-PS a été 
introduit, adapté et appliqué pour la solution optimale du problème d’écoulement de 
puissance dans un système électrique équipé de dispositifs SVC et TCSC. Ce problème de 
recherche des flux optimaux des puissances de circulation a été formulé comme 
d'optimisation combinatoire sous contraintes d'égalité et d'inégalité. L'algorithme proposé a 
été, ainsi, testé sur le modèle «test IEEE 30-bus » pour différents types de fonctions 
objectives. Les résultats montrent que l'utilisation combinée des dispositifs SVC et TCSC 
peut efficacement réduire le coût de production, améliorer la déviation de tension, 
diminuer les pertes actives et réactives et améliorer l'indice de stabilité de la tension dans le 
système électrique. Les résultats de la simulation montrent également l'efficacité de 
l’algorithme proposé pour la résolution du problème d’écoulement de puissance optimale 
pour un système électrique équipé de dispositifs FACTS. Les résultats des algorithmes 
proposés ont été comparés avec ceux d'autres méthodes récentes (citées en référence). 
L'évaluation comparative confirme, en termes de qualité de la solution, la robustesse de la 


méthode proposée pour tous les cas d'étude considérés. 
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5.3. Optimisation Multi objectif de l’écoulement de puissance par l'algorithme 


Multi-objectifs Jaya (MOJaya). 


Cette partie de travail présente les résultats d'optimisation de l'écoulement de 
puissance basée sur l'application de l’algorithme Multi-objectifs Jaya. Cinq fonctions 
objectives ont été optimisées individuellement et simultanément à savoir, le coût de 
combustible (cout), la déviation de tension (VD), l’amélioration de la stabilité de tension 
(Lmax), les pertes actives de transmission (Piosses), la sécurité du système (SD). Les cas 
étudiés dans cette partie de travail sont présentés dans le fableau 5.17. 

Le programme développé a été implémenté en utilisant le logiciel MATLAB. Les 
simulations ont été effectuées en utilisant l'approche MOJaya proposée avec n = 100, nA = 


100 et un nombre d’itérations maximal égale à 500. 


Tb.5.17. Résumé des cas étudiés. 


BCE ECO ES 


= JTE 


Mono-objective 
Mono-objective 
Cette étude aidera à identifier la gamme de variation de chaque fonction objective 


Mono-objective 


Mono-objective 
Multi-objective 
Multi-objective 
Multi-objective 
Multi-objective 


SUR 


5.3.1. Problème Mono-objectif. 


et aidera à analyser la meilleure solution de compromis plus tard. Les résultats obtenus de 
cette partie d’optimisation mono-objective sont présentés dans le tableau 5.18. De plus, la 
convergence des objectifs est représentée sur les figures 5.913. On peut remarquer sur ces 
figures que la fonction objective converge rapidement vers la valeur optimale alors que le 


terme de pénalité atteint zéro après quelques itérations. 
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Tb.5.18. Résultats optimaux trouvés pour les cas à objectif unique utilisant Jaya. 


Peu 
Pc13 


Va1 (p.u.) 


Vou 
Vos 


Ti (6-9) 


T2 (6-10) 


T5 4-12) 
T36 28-27 
QCio 
QC» 
QCis 
QC 
QC0 
QC 
QC: 
QC 
QC 
Cost ($/hr) 
VD (p.u.) 
L_max (p.u.) 
Ploss (MW) 
SI 


177.0386 
48.6817 
21.3210 
21.0966 
11.8740 
12.0001 

1.1000 
1.0807 
1.0541 
1.0619 
1.1000 
1.1000 
1.0216 
0.9000 
0.9645 
0.9530 
4.9998 
5.0000 
4.9955 
4.9999 
4.2670 
4.9998 
2.6187 
5.0000 
2.3064 


89.0808 
78.6206 
49.8306 
34.6289 
23.9941 
12.0077 
1.0248 
1.0143 
1.0127 
1.0071 
1.0441 
1.0004 
1.0646 
0.9010 
0.9574 
0.9699 
4.4080 
0.0000 
4.8290 
0.0773 
4.9988 
48611 
4.9784 
4.9206 
2.5858 


53.4323 
79.4068 
49.6937 
34.2490 
29.9464 
39.7724 
1.0991 
1.0925 
1.0868 
1.0783 
1.0998 
1.0999 
0.9791 
0.9063 
0.9746 
0.9437 
4.3023 
4.0689 
3.4300 
4.2433 
4.7195 
0.0426 
0.1934 
0.0277 
0.1018 


51.2433 
79.9994 
50.0000 
34.9999 
30.0000 
39.9999 
1.1000 
1.0932 
1.0752 
1.0822 
1.1000 
1.1000 
1.0526 
0.9000 
0.9836 
0.9686 
4.9936 
4.9955 
3.8780 
5.0000 
3.6541 
5.0000 
2.4078 
4.99388 
2.0461 


78.5505 
79.9934 
49.9996 
34.9847 
29.9922 
13.5290 
1.0999 
1.0914 
1.0584 
1.0704 
1.1000 
1.0948 
0.9799 
0.9797 
0.983815 
0.9794 
4.9830 
5.0000 
2.8708 
4.9847 
3.7408 
4.9775 
1.1458 
4.9916 
2.5062 


799.0343 | 907.2475 | 963.1270 | 967.0467 | 919.9918 


1.9737 

0.1260 
8.6121 
14.3861 


0.0935 
0.1488 
4.7626 


13.1574 
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2.0425 
0.1245 
3.1006 


13.0120 


2.0418 
0.1258 
2.8425 


12.6494 


1.7959 
0.1293 
3.6494 
11.2312 
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125 250 
Iterations 


x 10 


Cost 
— — — Penalty 


7 


Penalty 


375 500 


Fig.5.9. Fuel cost minimization. 
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Fig.5.10. Voltage deviation minimization. 
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Fig.5.12. Active power losses minimization. 
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Fig.5.13. SI minimization. 


À partir du tableau 5.18, on peut voir que la meilleure valeur du coût est de 
799,0343 $ / h et la plus faible est de 967,046 $ / h obtenue lorsque le Pisses St minimisé. 
La meilleure valeur de VD obtenue est de 0,0935 et la mauvaise valeur obtenue lorsque 
Lmax est minimisée. Pour Lx, la meilleure valeur obtenue est de 0.1245, alors que la 
mauvaise valeur est de 0.1488 obtenue quand VD est optimisé. En ce qui concerne les 
pertes, les meilleures et les mauvaises valeurs obtenues sont respectivement de 2,854 MW 
et de 8,6121 MW. Enfin, pour SL la meilleure valeur atteinte après optimisation est de 
11.2312 alors que la mauvaise (14.3861) est obtenue lorsque le coût est minimisé. 

Afin de montrer l'efficacité de l'algorithme de Jaya, les résultats obtenus pour le 
premier cas (puisque ce cas est le cas le plus étudié) sont comparés à de nombreux autres 
algorithmes comme le montre le tableau 5.19. On peut voir que Jaya est le meilleur parmi 


beaucoup d’algorithmes en termes de fonction objective trouvée. 
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Tb.5.19. Comparaison des résultats de simulation du CAS-1 avecdifférentes méthodes. 
799.0343 
799.0352 | Enhanced Colliding Bodies Optimization 
799.0353 | Improved Colliding Bodies Optimization 
799.0463 | Grenade Explosion Method 
799.0715 | Teaching-Learning-Based Optimization 


799.0760 | Backtracking Search Optimization Algorithm 


799.1116 | Biogeography-Based Optimization 
799.1821 | Improved Electromagnetism-Like Mechanism 
799.1974 | League Championship Algorithm 
DE 799.2891 | Differential Evolution 
SA 799.45 | Simulated Annealing 
PSO-PS 799.8723 | particle swarm optimization (PSO) and 
pattern search (PS) 
BHBO 799.9217 | Black-Hole-Based Optimization 
PSO 799.9738 | Particle Swarm Optimization 
EM 800.078 | Electromagnetism-Like Mechanism 
EADHDE | 800.1579 | Genetic Evolving Ant Direction HDE 
EADDE 800.2041 | Evolving Ant Direction Differential Evolution 
FPSO 800.72 | Fuzzy Particle Swarm Optimization 
IGA 800.805 | Improved Genetic Algorithms 
GAF 801.21 | Fuzzy Genetic Algorithm 
ICA 801.843 | Imperialist Competitive Algorithm 
EGA 802.06 | Enhanced Genetic Algorithm 
TS 802.2900 | Tabu Search 
802.376 | Modified Differential Evolution Algorithm 
802.465 | Improved Evolutionary Programming 
802.62 | Evolutionary Programming 
804.02 | Refined Genetic Algorithm 
804.853 | Gradient Method 
805.94 | Genetic Algorithm 


5.3.2. Problème Multi-objectif. 


Après résolution, dans la partie précédente, des cas d’écoulement optimal de 
puissance à objectif unique, on considère dans cette section les cas Multi objectifs. Puisque 


le coût est l'objectif le plus important, il convient de mentionner qu’il est considéré comme 
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une fonction objective dans tous les cas du traitement Multi objectifs. Les meilleures 
solutions de compromis trouvées pour les cas étudiés sont présentées dans le tableau 5.20. 
Ici, 1l est important de confirmer que toutes les valeurs des variables du vecteur de contrôle 
sont, pour tous les cas traités, dans leurs limites admissibles. Les ensembles optimaux de 
solutions trouvées de Pareto, correspondants à tous les cas considérés, sont illustrés sur les 


figures 5.14-17. 


Tb.5.20. Meilleures solutions de compromis « Multi Objectives - MOJaya ». 


174.7237 | 178.9573 | 129.0885 | 145.6301 

48.5100 | 48.5201 | 50.8800 | 43.9291 

18.3596 | 19.8210 | 29.7657 | 24.9058 

29.2951 | 18.6096 | 35.0000 | 35.0000 

10.0000 | 14.3773 | 30.0000 | 26.9082 

12.0000 | 12.0000 | 14.2854 | 14.3550 
1.0574 1.1000 1.1000 1.0380 
1.0385 1.0806 1.0854 1.0227 

0.9963 1.0594 1.0549 0.9650 

1.0000 1.0606 1.0701 0.9880 

Vu 1.0013 1.1000 1.0716 1.1000 
Ve 1.0360 1.1000 1.0706 1.0960 
Eee 1.0139 0.9978 1.0275 0.9000 
TER 0.9187 0.9000 0.9472 0.9000 
TS 1.0184 0.9686 0.9993 0.9048 
Pen 0.9783 0.9316 1.0220 0.9001 
QCio 3.9841 0.7551 3.4799 0.1773 
QC» 3.0661 3.7056 4.2978 0.0000 
QCis 5.0000 0.0000 0.9800 3.6568 
QC 5.0000 5.0000 5.0000 5.0000 
QC 5.0000 | 4.5824 5.0000 | 44096 
QCA 0.5857 5.0000 0.0000 5.0000 
QC: 4.9274 5.0000 2.0799 2.2283 
QC 3.8348 0.0000 5.0000 | 4.6058 
QC 4.5394 0.0000 33511 2.5925 

Cost ($/hr) 803.8962 | 799.7642 | 824.5752 | 816.1213 
VD (pu.) 0.1648 1.8657 1.1499 1.5126 
L_max (pu.) 0.1492 0.1254 0.1388 0.1314 
Ploss (MW) 9.4885 8.8853 5.6196 7.3283 
SI 14.3080 | 14.8357 | 12.8461 | 12.6369 
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Fig.5.14. Ensemble optimal de solutions de Pareto pour CAS-6. 
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Fig.5.15. Ensemble optimal de solutions de Pareto pour CAS-7. 
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Fig.5.16. Ensemble optimal de solutions de Pareto pour CAS-5. 
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Fig.5.17. Ensemble optimal de solutions de Pareto pour CAS-9. 
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Le CAS-6 donne la meilleure solution de compromis, qui peut être représentée par 
la paire (Coût, VD) égale à (803.8962, 0.1648). Cette valeur peut être comparée avec le 
CAS-I et le CAS-2 représentés par les solutions de valeur (799.0343, 1.9737) et 
(907.2475, 0,0935), respectivement. De plus, l'ensemble optimal de Pareto obtenu pour ce 
cas et illustré sur la figure 5.14 présente une bonne répartition entre les deux objectifs. Les 
points varient entre (800.7102, 0.7323) et (824.4804, 0.1551). 

Le coût minimal obtenu pour le CAS-7 est de 799.7642, avec une valeur de Lmax 
engendrée égale à 0.1254. Cette solution de compromis peut être comparée avec les 
valeurs obtenues dans le CAS-I et le CAS-3 représenté par les solutions de valeur (Coût = 
799.0343, Lmax = 0.1260) et (963.1270, 0.1245), respectivement. L'ensemble optimal de 
solutions de Pareto obtenu pour ce cas est illustré sur la figure 5.15. 

Pour le CAS-8, le coût total obtenu est de 824.5752 et les pertes actives engendrées 
égales à 5.6196. Cette solution de compromis peut être comparée avec le CAS-I et le 
CAS-4 représentés par les solutions de valeur (799.0343, 8.6121) et (967.0467, 2.8425), 
respectivement. La figure 5.16 représente les solutions de Pareto du coût et pertes actives 
optimaux pour ce cas. 

La meilleure solution de compromis trouvée pour le CAS-9 est de (Coût = 
816.1213 et SI = 12.6369). Cette valeur peut être comparée avec celles trouvées par le 
CAS-I et le CAS-5 représentées par les solutions de valeur (799.0343, 14.3861) et 
(919.9918, 11.2312), respectivement. L'ensemble optimal de solutions de Pareto obtenu 


pour ce cas est représenté sur la figure 5.17. 


5.3.3. Commentaire. 


Dans cette partie de travail, un nouvel algorithme multi-objectif Jaya a été 
développé afin de résoudre le problème d’écoulement optimal de puissance; lequel est 
formulé comme une optique d'optimisation combinatoire. L'algorithme développé a été 
implémenté et appliqué avec succès pour résoudre le problème posé sur le support test 
IEEE 30-bus; où neuf cas ont été étudiés (5 cas d’objectifs uniques et 4 cas multi- 
objectifs). 

Pour les cas à objectif unique, cinq fonctions objectives ont été considérées; à 
savoir, minimiser le coût total du combustible, améliorer le profil de tension, améliorer la 


stabilité de la tension, réduire les pertes de puissance active et améliorer la sécurité du 


117 


système. Pour les cas multi-objectifs, la minimisation du coût total du combustible est 
considérée avec chacune des autres objectifs. 

Les ensembles optimaux de Pareto obtenus pour ces cas ont été représentés pour 
montrer la performance de l'algorithme développé. Les résultats obtenus montrent la 


qualité des solutions lorsqu'une ou deux fonctions objectives sont considérées. 
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Conclusion. 


Le travail réalisé dans le cadre cette thèse est une contribution dans l’amélioration 
des processus de l’écoulement des puissances active et réactive dans les réseaux électriques 
en présence des dispositifs FACTS. Il consiste en un traitement combinatoire, orienté vers 
la recherche des relations optimales entre des variables de contrôle à même de minimiser 
une ou plusieurs fonctions « objectif » tout en satisfaisant un ensemble des contraintes 
d'égalité et d'inégalité. 

L'augmentation incessante des variables de contrôle impliquées, surtout en 
présence de dispositif FACTS, et la non linéarité des contraintes, rendent le problème 
d'OPF de plus en plus complexe; ce qui nécessite, pour atteindre les meilleures solutions, 
des méthodes performantes et robustes, comme celles très actuellement introduites et qui 
s’appuient sur des règles intelligentes dites meta heuristique. Dans ce travail, trois 
techniques méta-heuristiques sont utilisées pour résoudre le probleme d'OPF. Il s’agit de la 
méthode d’optimisation par essaims de particules (PSO), d’une méthode méta heuristique 
hybride qui combine l’algorithme des PSO et l’algorithme de PS et d’une nouvelle 
methode dite Multi-objectifs Jaya. 

L’application de ces méthodes, les analyses comparatives, les déductions objectives sur 

le plan des résultats des divers traitements de simulation montrent que : 

e Les méthodes metaheuristiques proposées peuvent contribuer à la résolution 
efficace des problèmes OPF. 

e Ces méthodes présentent une grande capacité de traitement, sont d’une grande 
robustesse et d’une très bonne flexibilité, 

e Les résultats obtenus sont très satisfaisants et très compétitifs par rapport à ceux 
obtenus par d’autres techniques, 

e L'utilisation des dispositifs FACTS améliore les résultats pour toutes les fonctions 


objectives traitées. 
En perspective; 


e Elargir l’application de ces méthodes à la résolution des problèmes d'OPF pour des 
réseaux électriques plus déployés (IEEES7 bus, IEEE 118 bus, IEEE 300 bus) en 
presence d’autres dispositifs FACTS; tel que SSSC, STATCOM, UPFC..., 
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Utilisation des méthodes metaheuristique pour déterminer les lieux optimaux 
d’emplacement des FACTS, leur type et la valeur de leur puissance; pour différents 
types de fonctions « objectif ». 

Appliquer ces methodes pour la résolution du problèmesd'OPF en tenant compte de 
la contrainte de securité ( Security constrained optimal power flow SCOPF j; 
laquelle contient des contraintes de sécurité supplémentaires qui assurent la sécurité 
du système non seulement en fonctionnement normal, mais aussi en cas d'urgence 
(panne d'une ligne de transmission et / ou transformateur ou générateur). Dans ce 
cas, on traite de deux principaux types de modèles SCOPF: l'approche préventive 


(PSCOPE) et l'approche corrective (CSCOPF)... 
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